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GEOLOGIE EN MINBOUW 


ON THE GRAVITY FIELD OF THE SUNDA REGION 
(WEST INDONESIA) 


B. J. COLLETTE 1 


SUMMARY 


Gravity observations at sea by Vening Mei- 
nesz and on land by the Bataafsche Petro- 
leum Maatschappiji in West Indonesia were 
subjected to a further correction for the influence 
of the sedimentary basins. This correction, the 
geologie reduction, has to be accompanied by a 
supplementary isostatic reduction. The data 
being scarce a 2-dimensional method, the so-called 
profile method, was applied in computing these 
effects. The resulting changes of the anomaly 
curves were small (fig. 4 to 8). 

The gravity field shows a mean positive 
anomaly of about 20 mgals. This regional anomaly 
was not taken into account.Important deviations 
from the picture of a negative zone, the Vening 
Meinesz zone, flanked on either side by positive 
zones — superposed on the regional field of 
+ 20 mgals — could be traced back to geologic 
causes (Part D. 

By means of applied mechanics (Part II) several 
solutions of the gravity field were tried. An 
upward force embodied by a root of light rock 
was imagined in an elastic crust floating in a 
fluid substratum. Only in case of very great 
thicknesses of the crust might sufficient agree- 
ment be obtained between the gravity field 
resulting from such a construction and the actual 
field. In addition, the strength of a crust of 20 
to 25 km thickness is not great enough to sustain 


INTRODUCTION 


The last one and a half century many 
investigations have been carried out on the 
geology and the geophysics of Indonesia. In 
the field of geology amongst others by the 
former Dienst voor de Mijnbouw (The Bureau 
of Mines) and by several commercial compa- 
nies which had and have their concessions in 
this country. 

With regard to geophysics, it was Professor 
Vening Meineszz who made many gravity 
determinations in submarines, while gravity 
surveys on land were carried out by the N.V. 
De Bataafsche Petroleum Maatschappij, and 
observations in other geophysical fields were 


1 Mineralogical-Geological Institute, University, 
Utrecht. 


the great stresses. From comparative studies by 
Vening Meinesz and Gunn it follows, however, 
that the crust’s thickness cannot be much 
more than about 50 km. So, if no other forces 
are presumed, the life time of the anomalies of 
the Vening Meinesz zone will entirely be deter- 
mined by the rate of reaction of the substra- 
tum. The half-life was computed from the data 
on the postglacial upheaval of Fennoscandia and 
ranged from 66 to 55.000 years. The latter figure 
will most probably not be enlarged when plastic 
instead of viscous flow is assumed. In order to be 
able to analyse the difference between the results 
of the computations on the viscosity of the sub- 
stratum by Van Bemmelen and Berlage (1934) 
— who assumed the viscous layer to have 
a finite thickness — and the other authors — 
who assumed the substratum to have an infinite 
depth —, the formula for the upheaval by 
Vening Meinesz (1937) was extended. Good 
agreement could be stated between the results 
obtained .with this extended formula and Van 
Bemmelen and Berlage’s. 

As the mentioned half life is much too short 
to account for the coexistence of negative zones 
in so many regions of the earth, the conclusion 
can hardly be avoided that considerable horizontal 
forces acting in the crust sustain the anomalies 
of West Indonesia. 


made by the Koninklijke Magnetisch en Mete- 
orologisch Observatorium (The Magnetical 
and Meteorological Survey) in Djakarta. 


To this considerable array of facts this 
paper will add no new data, except some 
gravimetric observations on land which were 
kindly put at my disposal for publication by 
the Bataafsche Petroleum Maatschappij, for 
which cooperation I here express my grateful 
thanks. 


In connection with these data several hypo- 
theses have been put forward on the genesis 
of Indonesia. On account of the special geologie 
position of this archipelago these hypotheses 
are far beyond the normal scope of a descri- 
bing regional geology and are therefore of 
the greatest importance for the understanding 


of tectonic processes in general. 
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The last few decades thorough-going dis- 
cussions have been held between the advocates 
of these hypotheses. In this study attempts 
will be made to offer some contribution to a 
few facets of this discussion. 

The first aim is to subject the gravity 
anomalies to a further correction. The usual 
gravimetric reductions assume a uniform 
density of the outer layers of the earth of 
2.67 or 2.6. This conscious idealisation naturally 
entails errors. It is desirable to apply an 
additional correction for this effect, for which 
correction the name geologic correction can 
be used. The appertaining isostatic reduction 
will be called the geologic isostatic reduction, 
the anomalies still remaining after this cor- 
rection, will be called geologie isostatic ano- 
malies. 

The methods by which these geological cor- 
rections are computed and its results are 
discussed in Part I. 

Secondly, it will be tried by means of 
applied mechanics to derive some limits 
within which hypotheses on mountainbuil- 
ding have to be confined. 

The purpose of the geophysical science 
with regard to the solid earth can be defined 
as describing the physical processes in the 
earth and determining the physical proper- 
ties of the material of which the earth has 
been built up and their spatial distribution. 
This is attempted with the aid of seismology, 
gravimetry, magnetometry and electrometry 
and by means of physical considerations of 
geologice and geomorphologic data. 

In these ‘purposes geophysics and geology 
can gain much by cooperation; scientists of 
both sides have the duty to try to remove 
contradictions between the results of these 
two branches of earth’s science, 

As to geophysics, this task begins with a 
clear statement of the consequences for 
geology of the geophysical data. 

In Part II we find some of these confining 
computations and statements that are based 
on the results obtained in Part I. 


PART I 
General remarks 


In order to be able to compare the various 
gravity data mutually and with respect to 
the theoretical value of gravity for the area 
in question, it is necessary to apply some re- 
ductions. 

First of all there is the free-air reduction 
to eliminate the effect of the different eleva- 
tions of the stations. 


Next, the topographic reduction is applied, 
in which a correction is .made for the irre- 
gular effect of the topography 2. It is assumed 
that the densityo of the topographic masses 
has a constant value, namely 2.67 or some- 
times 2.6, 

The isostatic reduction is the third. On the 
ground of deflections of the vertical, Airy 
presumed an isostatic compensation of the 
topography caused by a floating of the earth’s 
cerust in a substratum. This implies that in 
case of (isostatic) equilibrium at a certain 
depth in the earth the pressure is constant 
at that level. 


This conception has been confirmed by 
many gravimetric observations. Several me- 
thods have been developed for computing the 
effect of such a compensation. The first me- 
thod, still used by the U.S. Coast and Geo- 
detic Survey, is: 

a. the method of Pratt-Hayford. To obtain 
a simple formula, Hayford based his tables 
for the effect of the isostatic compensation 
on Pratt’s hypothesis. In this hypothesis the 
topographic relief is assumed to have 
originated from expansion of vertical columns 
of equal height. So the compensation is 


2 This reduction is occasionally identified with 
the Bouguer correction. The Bouguer correction, 
however, is only a first approximation of the 
topographic correction, with the assumption that 
the effect of topography can be equalized to that 
of an infinite plane slab, having a thickness 
equal to the elevation of the station. For this 
we have the formula: Ag = 2,rkhso, 


in which Ag = the attraction of this plane slab 
k = Newton’s constant 
and h = the elevation of the station 


or Ag= 0.1118 AR mgals, in which} is expressed 
Ss Ss 
in meters. 


By a mere coincidence this formula is in 
accord with the first and most important term 
of the formula for the topographie correction 
with Bullard’s method, in which the attraction 
of a cap with radius of 166.7 km and of a 
certain thickness is computed. Thus is obtained 
the formula. Ag = 0.1118 h, + B mgals, in which 


B is a slight additional correcetion dependant on 
the elevation of the station. 


The deviations of actual topography from a 
cap is computed separately. This latter correction 
is always positive, since inward as well as outward 
bulgings have a diminishing effect on the gravity 
value at the point of observation (see Lejay, 
1947, p. 65 f£.). Regardless of the correction 
method followed the name Bouguer anomaly is 
used for the remaining anomaly. The effect of 
topography of the regions outside this cap is 
computed together with the isostatic effect. 


thought to be equally ‘distributed over this 
column. 


b. Heiskanen later on developed tables 
for Airy’s idea, that is, a floating crust on 
a substratum behaving like a fluid. For his 
tables he assumed several values for the 
crust’s thickness. In Airy’s conception the 
crust can be most easily imagined as icefloes 
of various thicknesses, floating in water, If 
we assume that these icefloes form one whole, 
this way of representing leads automatically 
to the conception of regional compensation. 
For not every irregularity of the upper sur- 
face of the floating ice is followed by its 
replica on the lower surface with an ampli- 
tude corresponding to the ratio of the densi- 
ties of-ice and water. Stresses may therefore 
occur in the floating body that are conside- 
rably greater than with strictly local isostatic 
compensation. 


c. For this expansion Vening Meinesz 
developed tables. His object was to obtain 
in this way material for comparison about 
the degree of regionality, that is, to ascertain 
whether the compensation has a strictly local 
character, or whether it is spread out hori- 
zontally beneath the topographic element in 
question and to what degree. His tables are 
worked out with various degrees of regiona- 
lity and with various values of the thickness 
of the crust. 


The geological reductions 


It is widely known that after all this treat- 
ment often very considerable anomalies re- 
main. This is partly due to the circumstance 
that the isostatic compensation is not every- 
where perfectly complete; this again can be 
caused by the presence of forces in the earth 
that give rise to deviations of the isostatic 
'equilibrium or interfere with its readjust- 
ment. Another cause must be looked for in 
inhomogeneities in the outer layers of the 
earth. 

In order to arrive at a more quantitative 
analysis of the remaining anomalies, compu- 
tations can be made about the influence of 
known geologic structures and their supple- 
mentary isostatic compensation. 

We want to stress that one had better 
start from the isostatic anomalies. Vreugde 
(1935) and Hofman and Van Bemmelen 
(1952), for instance, based their consider- 
ations on the gravity field — i.c. of North 
Java and the Western Mediterranean — on 
the Bouguer values. However, isostasy is a 
solid fact. Now, if a minor irregularity in the 
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gravity field of a flat area is to be investiga- 
ted, one might start from the free-air ano- 
malies. In such a case these anomalies are 
to a fairly high degree in accord with the 
isostatic ones, 

But when larger changes in horizontal 
sense of gravity are analysed, like in the 
cited publications, and the Bouguer values 
are used, then the effect of the isostatic com- 
pensation of the area itself and of the sur- 
rounding regions on this horizontal change, 
the horizontal gradient of gravity, is neglec- 
ted. In this way incorrect conclusions might 
be reached and so in principle this method 
has to be rejected. 


It is true that with the isostatic reduction 
it is often not exactly known which system 
of correction is the right one in the case in 
question. This very uncertainty may be a 
warning for too detailed and unjustified pro- 
nouncements. 

In his computations on volcanic islands 
Vening Meinesz (194la) was one of the 
first to introduce the effect of geologic struc- 
tures. Vening Meinesz started from the 
fact that the density of the rock of volcanic 
islands is greater than the theoretic value 
2.67 and involved this in his computations. 
Further, Evans and Crompton (1946) were 
the first who systematically computed the 
effect of sedimentary basins over a large 
area. They did so for a part of Burma and 
India and it was cliearly demonstrated how 
great the effect of sedimentary deposits of 
these basins and their compensation can be. 

So far such computations have not been 
performed for Indonesia. In the following 
pages this gap with regard to the Sunda 
region is filled. 

For this purpose a geologic and the apper- 
taining geologie isostatic correction is com- 
puted for some observations at sea of 
Vening Meinesz (1934 and 1948) and obser- 
vations on land by the B.P.M., who kindly 
released some data for this publication. 

These observations on land are represented 
in Table I. They are supplemented by obser- 
vations of the same company in Java, pu- 
blished by Vreugde (1935), who in his turn 
made a supplementary use of Lejay’s ob- 
servations (republished in 1936). 


The gravimetric data as well as the geolo- 
gic ones are not suitable to such a widely 
ranged plan as Evans and Crompton's. 
A 2-dimensional way of approach, the so- 
called profile method, will therefore be ap- 
plied. This method is less laborious than the 


274 


TABLE I — GRAVITY DATA ON LAND OF THE B.P.M. 


—————————————————————————RRREE 
Elevation Observed Bouguer 
No. Latitude Longitude Ber g anomalies 
a ET ken ar il nern lin en nenn. 
Sumatra 0= Zu 
3°17.0’N 99°14 3’E 51 978.054 — 

: 02.3 40 .6 12 .082 + 21 

9 2 20.0 01.3 1072 977.809 az. 

10 1 46.1 98 48.2 323 .980 —0 

15 001.9 100 44 .6 149 973.049 + 30 

16 18 .7S 222 922 977.849 u 

17 37.6 07.0 2 978.063 + 14 

18 1 20 .3N 102 09 .3 3 .068 + 77 

19 0 32.3 101 26 .4 6 .063 + 15 

20 57 .6S 100 36 .6 955 977.855 — 8 

21 40 .1N 103 06 .4 2 978.066 + 17 

22 09.9 102 33 .8 3 .075 + 26 

23 29 .7S 04 .4 18 .055 + 9 

24 1 05.3 101 53 .9 90 .064 + 30 

26 1 05.5 104 11.7 2 .078 + 28 

27 1 42.3 103 06 ..3 25 .079 + 30 

28 241525 01.3 78 .062 + 20 

29 1 54.4 104 21.0 2 .106 + 51 

30 2 45.0 103 33 .8 12 .074 + 16 

31 322.1 03 .4 58 .090 + 34 

32 347.2 102 19.7 17 17 + 9 

34 2 25.5 105 332% 2 .092 + 35 

35 2 48.1 104 44 .8 3 an + 50 

36 3 38.9 34 .6 20 117 + 51 

37 4 04.1 01.1 92 .091 + 33 

Java v0o= 260 

4 6°11.0’S 106°50 .1’E 6 978.159 + 51 

5 37.1 48 .3 280* 111 + 50 

6 59.5 32.6 15% 250 + 128 

7 714.9 25.8 ® _ + 225 

19 6 45.0 IIE0TES 12 127 + 9 

20 7 34.3 110 49.8 90 .097 — 23 

21 744.8 52.7 105 .195 + 73 

22 8 01.2 56 .8 160 .213 + % 

23 8 14.2 111 03.6 196 222 + 106 
a a Fre EURE Dr EEE Fe EA En a nn ion Elson... 
Column 1 All numbers represent Holweck-pendulum-stations. 

„ 2/3/4 The observers plotted the stations on the excellent 1:50.000 or 1: 100.000 topo- 
graphic maps with contour lines, issued by Top. Dienst, Batavia, and the given 
values should therefore be accurate. The values marked with an asterisk are not 
the original values from the observers, but have been taken from prints in the 
office, since most of the records of the Java-surveys are missing. They are 
therefore somewhat less accurate. 

„ 3 Java: 

The g-values have been calculated back from the Bouguer-anomaly figures, for 
the same reason. The survey depends on the bases Batavia eg = 978159 and 
Bandung 978,988. 
Sumatra: 
The g-values depend upon: 

De Bilt 978.267.5 

Batavia M.I.. 978.159 

Bandung 978.988 

6 


gas mus anomalies have been computed, using Cassini’s 1930 International 
ormula. 


3-dimensional one büt naturally also less 
accurate. In areas of a pronounced linear 
character and for a limited range, it may be 
a good approximation, however. 

Figure 1 shows the position of the profiles, 
namely three over Sumatra and two over 
Java, as well as of the used observations. 


The isostatic reduction of the gravity data 


The observations on land by the Bataafsche 
Petroleum Maatschappij have been supplied 
as Bouguer anomaliess. By means of the 
”Isostatic Isocorrection Line Maps”, a publi- 
cation by the E.A.E.G. (see list of references), 
it was possible to reduce the anomalies iso- 
statically; however, only with the system 
worked out in this publication, for a crust’s 
thickness T = 25 km, and a regionality R = 
160 km. 

In our computations we must confine our- 
selves to these regional isostatic anomalies ®. 
In order to have some comparative material 
available in spite of these limitations, we have 
computed for one of the profiles (profile IV) 
also the geologic local-isostatic anomalies. 

Basing our computations on these regional 
isostatic anomalies, it is useful to realise 
where we can expect that these anomalies 
are a correct or less correct apnroximation. 

In the first place there is the problem of 
the transition from continent to ocean, for 
which Vening Meinesz (1948) demonstra- 
ted that this transition in the normal case is 
locally compensated. So here we may expect 
a considerable deviation. 

For the Vening Meinesz zone on the 
other hand, it is quite probable that a regional 
form of reduction is a good approximation. 
This zone has been forced out of equilibrium, 
so there are important stresses in the 
crust and since we must assume that some 
time is necessary to make these disappear, re- 
gional reduction based on physical deforma- 
tion laws will better fit in with this picture. 

As for the geological older areas, of those 
it must be assumed that the compensation is 


3 When comparing the values stated here for 
Vening Meinesz’s stations at sea (1934) it will 
be seen that there are some little differences. 
This is due to the Browne correction applied 
later on, a correction for deviations of the second 
order caused by accelerations of the ship in 
which the apparatus is adjusted. For this effect 
Vening Meinesz (1953) developed the com- 
plete formulas. Those corrections are listed in 
a later publication by Vening Meinesz (1941b, 
1948). The anomalies thus corrected were also 
used by De Bruyn (1951) in compiling his 
gravimetric map of SE Asia. 


275 


of local character, Fortunately, these areas — 
e.g. the Java Sea — are for the greater part 
fairly flat, so that it makes little difference 
what system is made use of for the reduction. 


The position of the profiles 

The position of the profiles was chosen 
according to the location of Vening Mei- 
nesz’s stations. They nearly coincide with 
his profiles 2, 3, 4, 6 and 7, on the understan- 
ding that our profiles form straight lines. 
Should comparison be made, this should be 
taken into account. The profiles are correspon- 
dingly numbered with Roman numerals II, 
III, IV, VI and VII. 


The sedimentary basins 


In the first instance some computations are 
given for the influence of the sedimentary 
basins on the gravity field. The data on the 
form of the basins and their extension are 
still rather scarce. No doubt, much has be- 
come known about the structure of the basins 
through the deep-borings of the oil-com- 
panies, but little has been published. For this 
reason we have often been compelled to take 
two depths for the basins in our computa- 
tions; in due course a choice may be made 
between them. 

When drawing the geologic profiles, we 
have mainly based ourselves on Van Bem- 
melen’s standardwork ”The Geology of In- 
donesia” (1949) and also on personal com- 
munications of geologists who have been 
working in Indonesia. Another source, chiefly 
important for the form of the basins was an 
article by Dufour (1951). 

It should be emphasized that the geologic 
profiles do not pretend to represent the cor- 
rect shape and extension of the sedimentary 
basins. On the contrary, it is fairly certain 
that locally the basin floor is disturbed tec- 
tonically, and moreover, it is doubtful whether 
it had such a regular curved shape from the 
beginning. The effect of certain disturbances, 
however, will be discussed when the remaining 
anomalies are considered, since actual data 
on this point are poor and scarce. 

With regard to the above, it is therefore 
of little importance to enter into all details 
of the geological profiles. Still it should be 
pointed out that the thickness of the. sedi- 
mentary layer of the outer arc naturally is 
but a speculation. On the island of Nias the 
total tertiary sedimentary column is of about 
4000-5000 meters thickness (Van Bemme- 
len, 1949, p. 164 ff.). The structure of the 
island is roughly an anticline, therefore the 
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sediments in the centre have not got this total 
thickness. Near Benkulen, the upper Pliocene 
is exposed (Van Bemmelen, 1949, p. 115) with 
an estimated thickness of 800—1000 meters. 
These are the only data and for this reason 
two thicknesses have consistently been taken 
for the sediments of the outer arc. 


For the ocean floor proper a sedimentary 
layer of 500 meters was assumed. The gravi- 
metric effect of such a layer and its compen- 
sation is small of course. 


The density of the sedimentary deposits 


In this paragraph the density of the sedi- 
mentary deposits of the basins is discussed. 
The only densities that are known of the 
Sunda region were published by Willems 
(1940). Of the more than 1200 analysed rock 
specimens only of less than a quarter has the 
density been determined. And among them 
there are not more than ten specimens of pure 
sedimentary nature. 


With regard to the oil area in Java, 
Vreugde (1935) starts from a density contrast 
of 0.5. However, this figure is not based on 
direct observation. We are therefore obliged 
to look for data of other regions. 


Our main source of informaticn is then the 
paper by Evans and Crompton (1946). For 
Recent to Tertiary rock these authors assumed 
densities of from 2.0—2.67 above 5000 feet and 
of from 2.1—2.67 if below it. 


In a detailed gravimetric survey in England 
Cook et al, (1952) arrived at values of from 
2.4—2.7 for the density of Triassic to Cambrian 
rocks. 


Hofman and Van Bemmelen (1952) neglected 
the influence of sedimentary basins, which in 
fact are not very extensive in the profile dealt 
with by them. 

Reich and Zwerger (1943) give values 
for the density of from 1.74—2.69 for sedimen- 
tary rocks of different character of pliocene 
to carboniferous age in Germany. 

Ewing in a paper not yet published * con- 
cludes to a value of 2.3 for the density for 
still unconsolidated sediments in and near the 
Porto Rico trough. 

From all these data we have decided to 
carry out computations with two densities for 
the sedimentary basins in the Sunda region, 
namely 2.3 and 2.5, while, moreover, we 


4 This paper has been published in the February 
issue of the Bull. Geol. Soc. America, Vol. 65, 1954. 
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combined the first value with the maximal 
depth and the second value with the minimal 
depth of the basins. 


The so-called profile method 


In dealing with a gravimetric problem with 
the socalled profile method, we started from 
the formula for the vertical component of the 
attraction of an infinite, homogeneous, hori- 
zontal elementary prism 


dx dz 
er) 

dgz x°+ z? 

(1% 
in which 
we the vertical component of 
the gravity 

k = Newton’s constant 

x = the horizontal direction per- 


pendicular to the prism 


z —= the vertical and 


dx dz = dS = the cross-section of the prism 
(Lejay, 1947, p. 205). 


When x and z in this formula are sub- 
stituted by Randp, 


in wich 
x—= R sin p 
z=Rcosp 

and dS is substituted by rdpdr, we get 
dg, = — 2 ko dz dp 


(2). 


On the ground of this formula it is possible 
to construct a diagram in which the attraction 
of each compartment is the same if dz and dp 
are taken as constant values (fig. 2). By 
counting the separate compartments and parts 
of them by wich 'the diagram covers the 
profile part to be computed, the total at- 
traction of each cylinder with arbitrary section 
can be obtained. 

On account of the problems in drawing, 
the vertical scale of the geologic profiles had 
to be exaggerated five times, and so the 
diagram had to be transformed into this scale 
too. This entails no changes in the principle. 


The line of density contrast 2.67/3.27 


Since we start from regional isostatic re- 
duction with a regionality R = 160 km, it is 
further necessary to find a method for draw- 
ing in profile the line of the supposed density 
contrast 2.67/3.27 that follows from this scheme 
of reduction. (In case of local isostatic 
compensation this line is simply at a depth 
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2.67/(3.27—2.67) —= 4.45 times the topographic 
height or (2.67—1.03)/ (3.27—2.67) = 2.73 times 
the sea-depth.) 


Fig. 2 — Diagram for computing the gravimetrical 
effect of 2-dimensional structures. 


On page 69 Vening Meinesz (1946) pre- 
sented in a table the subsidence expressed 
in meters caused by a topographic elevation 
of 1000 meters for the 2-dimensional case, with 
values for R of 58.1, 116.2, 174.3 and 232.4 km. 
In computing this subsidence it was not 
necessary to simplify the curve of downward 
bending, as had to be done in the 3-dimensional 
case (Vening Meinesz, 1940, p. 12). So here 
also negative values appear at the outer 
boundaries, i.e, bulgings. 

The degrees of regionality are not directly 
suitable for us; moreover, the anomalies used 
by us are computed with the aid of the 
simplified 3-dimensional method. Therefore 
attempts were made to find another method, 
with as point of departure the curve of down- 
ward bending on page 13 of the_cited 
publication (1940), that served as a basis for 
the 3-dimensional tables. This curve was in- 
tegrated in order to arrive at the corresponding 
curve for the 2-dimensional case. Next the 
surface enclosed by the curve thus obtained 
and the horizontal, was divided into 31 columns 
of equal bases (fig. 3). The total surface of these 
columns was put equal to 4.45 respectively 
to 2.73 times the cross-section of an infinite 
bar, 100 meters high, 10 km wide, out of 
which the values for the ordinates of the 
points of the curve could be found (see 
Table II). 


With the aid of these figures a table was 
composed, through which by means of simple 
counting, going along the profile, the depth of 


rn 


Cross-secf/on of 
a infinite bar 


Verfies! scale 
1250x ezagerrsfed 


Ile 


Fig. 3 — The regional isostatie compensation of 
an infinite bar. 
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the line of density contrast can be computed 
at any desirable point. This table can also 
be applied in computing the compensation of 
geologice structures with deviating densities, 
by transforming the structure into a pseudo 
topography with the aid of the formula 


h, = (1 - 09/o,) .h 


(3), 
in which 
h, — the pseudo height 
ne the density of the sur- 
z rounding rock 
0 = the density of the structure 
to be corrected 
and A = its vertical dimension. 


The geologic isostatic anomalies 


By means of the methods and assumptions 
discussed in the preceding paragraphs, it was 
possible to compute the geological corrections. 
The results are presented in Table III. For the 
drawing of the curves (fig. 4-8) the mean 
values were taken of the geologie isostatic 
anomalies belonging to the assumption of 
maximal and those belonging to the assumption 
of minimal sedimentary cover. 


In drawing the curves of the normal isostatic 
anomalies supplementary use has been made 


TABLE II — RELATION BETWEEN TOPOGRAPHY AND COMPENSATION IF R — 160 KM 
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TABLE III — THE GEOLOÖGIC ISOSTATIC ANOMALIES. 
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TABLE III — Continued. 
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*) The italicized numbers indicate the stations of 
Vening Meinesz. The normal ones are B.P.M.-obser- 
vations (see table I). 


of the isogam map by De Bruyn (1951) and 
of the Bouguer anomaly map of N Java by 
Vreugde (1935). The values derived from 
De Bruyn are indicated in the curves by 
the letter B, those derived* from Vreugde by 
the letters A to G. 

Some stations are too far removed from the 
profile line to serve directly as data for the 
curve in question. In that case a cross sign 
(x) indicates the isostatic anomaly belonging 
to that point. 


Examination of the results 
An examination of the curves tells us that 
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gives the 
geologic isotatic anomaly pertaining to the assumption 
of a thin sedimentary cover with density 2.5, the 
lower part gives the anomaly if the sedimentary 
cover is thicker with density 2.3 (see text). 


**) The upper part of the sign or 


the correction for the sedimentary basins 
entails only slight modifications of the total 
picture. Moreover, the largest changes are 
not only determined by the geologic reduction 
itself, but also by the reduction method applied 
for the supplementary isostatic compensation 
of the basins. If we had chosen a local reduc- 
tion, the changes would have been smaller. 
This is clearly illustrated in profile IV, where 
it is seen that the geologic local isostatie 
anomaly as regards the continental part at 
least, virtually coincides with the norma 
regional isostatic anomaly. 


For the rest we may conclude that the 
introduction of a geologie correction for the 
sedimentary basins does not throw new light 
on the whole of the gravity picture and entails 
no essential changes for the picture — as fixed 
by Vening Meinesz — of a negative zone 
with positive zones on either side. 


Further remarks on the geologic isostatic 
anomaly curves 


We shall now try to give a more detailed 
description, of the remaining anomaly picture 
on the basis of more general considerations. 


In the first place it should be remarked that 
the entire level of the curves is too high, that 
is, the positive anomalies are not compen- 
sated by the negative ones. Vening Meinesz 
(1934, p. 109) computed that the average 
anomaly amounted to + 20 mgals for all his 
observations at sea in Indonesia. This figure 
will not change considerably if the observa- 
tions on land are involved in such a compu- 
tation. 

So in Indonesia we have to do with a large 
regional positive field. It will not be tried 
to explain this. Similar regional fields are 
met with in India and in the North Atlantic. 
Thus, assuming this positive field to be of 
regional character, in accounting for the 
anomalies we shall have to start from an 
auxiliary level at +20 mgals, whereby it 
should be borne in mind that this is naturally 
but an approximation. 

Secondly it can be seen that there are still 
important deviations from the picture of a 
negative zone with positive zones on either 
side: 

1. In profile II there is a negative near the 
Toba Lake, already very important in the 
curve of normal isostatic anomalies, still more 
important in that of the anomalies corrected 
for the influence of the sedimentary basins, 
the geologic isostatic ones, and in absolute 
value even greater than the negative of the 
Vening Meinesz zone. Also in profile III 
a similar negative occurs, and to a less extent 
in profile IV, in which the minimum is only 
relative. 

%. In the profiles II, III, IV and VI, there 
are deviations of a shorter wave length on 
the land part. 

3. Profile VI shows a very high positive 
north of the Vening Meinesz zone. 

4. Profile VII gives an important negative 
near the Kendeng zone, which in itself consists 
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of two negative parts divided by a positive. 
ad 1. The region of this (relative) minimum 
is known in the geologic literature as the 
Semangko zone, a graben that stretches from 
Sabang to the north of Sumatra, to the bay 
of Semangko in the extreme south of Suma- 
tra (Van Bemmelen, 1949). Van Bemmelen 
believes the Semangko zone to be a graben on 
the top of a geanticline. Vening Meinesz 
considers it to be a shearzone resulting from 
the assumed movement from E Sumatra to 
the SE, which assumption is supported by the 
displacements after an earthquake observed 
by Muller (1895). In profile II the minimum 
shows two peaks; the SW peak is somewhat 
more important than the NE one. Also De 
Bruyn’s map (1951) clearly demonstrates 
that this minimum actually consists of two 
confluent negatives (see fig. 9). The Toba 
Lake depression (Van Bemmelen, 1949) 
could be the cause of the NE negative. In the 
sub-soil light rock or partly emptied magma 
chambers may still be present. 


EURE (-2 2) m gals 
777 From O fo - 2dmgals 
5 +25 fo Omgsls 


Y+25mgals 


Fig 9 — Regional isostatic anomalies (R = 160 
km, T = 25 km) of the Toba Lake Region (after 
De Bruyn, 1951). 


The assumption that the Semangko zone 
is a topgraben and the assumption that it 
represents a shearzone, both offer sufficient 
possibilities to account for the minimum. An 
attempt to come to a more quantitative 
estimation of the causes of the negative would 
be idle, since the number of solutions is 
geometrically infinite. 

ad 2. The deviations of a smaller wave 
lenght, also occurring like dome-shaped 
maxima and minima on De Bruyn’s map 
may be explained by assuming inhomogeneities 
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in the crystalline basement complex or in the 
sedimentary basins themselves, possibly 
accompanied by irregularities in the basin 
floor. An indication for the correctness of this 
suggestion may be found in the structure of 
the island of Banka (gravimetrically a relative 
minimum), where the tingranite is possibly 
surrounded by denser rock like the rock near 
Muntok (Western Banka) indicates (Wester- 
veld, 1939). 


ad 3. The considerable positive near the 
Wijnkoopsbaai (north of the Vening Meinesz 
zone in profile VI), greater than 150 
mgals, cannot be explained subcrustally on 
account of its large gradient. This positive is 
superimposed on the positive zone next to the 
Vening Meinesz zone. The positive seems 
to have a zonal character and may be traced 
all along the coast of South Java (see profile 
VII) and partly along Sumatra (see profile 
IV). The Wijnkoopsbaai zone may further 
appreciably be seen continuing in profile II 
and III, where the gradient is much smaller, 
however. 


A geologie indication for the explanation 
is to be found in the presence of gabbro, 
peridotite and serpentine in the Pretertiary 
of the Tjiletuhbay, a part of the Wijnkoops- 
baai. 


An analogy of the Wijnkoopsbaai zone is 
found in the Ivrea zone south of the Alps, 
described by Niggli (1946) and Vecchia 
(1952). Here too, that is, in a tectonically 
corresponding place, highly positive anomalies 
are found combined with exposures of denser 
rock. Niggli carried out some computations 
with more or less satisfactory results, 
Vecchia only presumed that the positive 
anomalies of the Ivrea zone are caused by a 
large dyke of dense rock. This dyke could 
have originated from the substratum during 
the orogenesis, 


Another explanation might be that we have 
to do with the edge of an isostatically com- 
pensated area with denser rock, formed 
during a previous stage of development. 
Gravimetrically the former explanation repre- 
sents a heavy vertical plate—which will 
be regionally compensated — the latter a 
horizontal one, and if the density surplus 
would be confined to the outer edge or upper 
edge respectively of these plates, the solution 
represents a horizontal cylinder. All these 
explanations are illustrated in fig. 10a, b and c, 
in which also the pertaining isostatic compen- 
sation is schematically indicated. 


Ad 4. On De Bruyn’s map the negative 
of the Kendeng zone stretches from Street 
Madura to south of Semarang. Further to the 
West, into West Java, a less clear continuation 
of the zone can be traced. In the plain of 
Djakarta, however, the anomalies are again 
positive (profile VI), just as over the idiogeo- 
syncline of East Sumatra (profiles II, III and 
IV), i.e. also positive with respect to the 
auxiliary level on +20 mgals. Thus a remark- 
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Fig. 10 — Possible causes of the positive of the 
Wijnkoopsbaai zone: 


a — vertical dyke, regionally compensated; 


b — horizontal layer, locally or regionally 
compensated; 
ce — horizontal cylinder, locally or regionally 
compensated. 


able difference appears between the idiogeo- 
syncline of East Sumatra and the plain of 
Djakarta on the one hand, and the Kendeng 
zone in East Java on the other. As to the 
geology of these areas it can be remarked that 
the Kendeng zone is situated S of a secondary 
range of volcanoes (for instance the Gunung 
Muriah and the Gunung Lasem). Van Bem- 
melen (1949, p, 730 and fig. 378; see also 
1954, fig. 51) describes this volcanic range as 
a part of the Pulu Laut system, that has to 
be considered as a southward extension of the 
Philippine system and that overlaps the older 
trendlines of the Sunda orogenese. For the 


rest there are no important differences in the 
geology of East Sumatra and North Java; for 
instance, the depth of the idiosynclines of both 
areas does not differ much (Umbgrove, 
1949, p. 44), The gravimetric difference now 
poses the following problem. If geosynclinal 
subsidence is a result of compression, we must 
expect negative anomalies above a geosyncline, 
and in that case the Kendeng zone is in normal 
condition. Assuming the depth of the idiogeo- 
syncline at 5000 meters (Umbgrove, 1949), 
this would yield an anomaly of more than 
—100 mgals, or, in consequence of the regional 
positive field, an anomaly of about —80 mgals. 
In East Java the anomaly is actually of that 
order of magnitude. If the root beneath Su- 
matra is assumed to be still present, some 
cause must eliminate the influence of the root 
and make the anomaly in Sumatra -20 
instead of —80 mgals. Such a cause is difficult 
to imagine and the conclusion seems to be 
justified that physical or physico-chemical 
processes made the root beneath the plain of 
Djakarta and the idiogeosyncline of East Su- 
matra disappear. Again assuming that the 
idiogeosynclines have a depth of about 5000 
meters and assuming a density contrast be- 
tween the crust and the substratum of 0.6, 
then this root could have formed mountains 
of about 1 km high. Now that neither these 
mountains and, considering the gravimetric 
anomalies, nor this root is present, the con- 
clusion may be made that the assumed pro- 
cesses of physical or physico-chemical nature 
are so important, that they may prevent 
mountain building. 


If, however, geosynclinal subsidence has to 
be seen as a result of changes in the sub-soil 
rock, by which the crust in situ becomes denser 
(positive anomalies) and subsides, the 
condition of East Sumatra and the plain of 
Djakarta would be normal (again with respect 
to the regional positive field), i.e. a slightly 
positive field could be expected over the 
geosyncline, and the Kendeng zone would 
represent a different case. 

A solution of the problem will not be sought 
here but it may be looked for in connecting 
the gravimetric data with the presence of the 
above-mentioned volcanic range. 


We need not be surprised that the negative 
of the Kendeng zone in itself is again divided 
in two parts, if we consider the geology of 
the zone. Beneath the dividing relative maxi- 
mum there is the Kendeng anticline, seen 
from tectonic angle un upfold of consi- 
derable dimensions. The data at our disposal 
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are not sufficiently detailed to undertake an 
attempt to a quantitative analysis starting 
from the geology, but Vreugde (1935) succeeded 
in giving a (geometric) solution of the gravity 
field, to which may be referred. 


The general pattern of the gravity field 
of the Sunda region 


Although some question-marks remained in 
the preceding paragraph, the greater part of 
the irregularities could be traced back to the 
geology. We are therefore justified in putting 
the gravimetric problem of West Indonesia as 
one of a negative zone flanked at either side 
by a positive zone, even though we do not 
know exactly in how far the positieve anoma- 
lies are to be regarded as being part of the 
vast positive field and of the positive zones 
sensu stricto. 

The picture positive-negative-positive dis- 
tinctly and symmetrically dominates in all 
the gravimetric profiles, be it that the negative 
anomalies in the profiles over Sumatra are 
weaker than in the profiles over Java. (It 
should be remarked that off Sabang the 
Vening Meinesz zone again aggravates, as 
appears from Vening Meinesz’s profiles 
and from De Bruyn’s map. A connection 
with the negative of Burma via the Andamans 
(—76 mgals) and the Nicobars seems likely.) 

In Part II the question will be dealt with 
in what way such an anomaly picture can be 
explained, in which the accent will be put 
on the question of the durability of the phy- 
sical configurations to be sketched. In many 
regions of the earth zones of negative anoma- 
lies are found, in SE Asia, the Antilles and 
in the Mediterranean. The assumption seems 
justified to consider these zones of negative 
anomalies as a stage in the orogenesis prior 
to the upheaval. From the geologic data it 
immediately follows that upheaval occurs only 
once in many millions of years in the same 
region. 

Supposing now the lifetime of a negative 
zone such as the Vening Meinesz zone to 
be only some tens of thousands of years, the 
coincidence of the presence of negative zones 
in so may regions would exceed the limits of 
probability, unless a reason could be found 
for such a world-wide occurrence just now 
in the holocene. 

Since such a reason is difficult to assume, 
we are obliged to conclude that the lifetime of 
{he anomalies of negative zones is geologically 
important, that is, anyhow considerably larger 
than for instance the order of magnitude of 
50.000 years. 


284 
PARIS 


The mass deficiency of the Vening 
Meinesz zone 

The cause of the negative anomalies of the 
Vening Meinesz zone must naturally be 
looked for in a deficieney of mass. This 
deficieney may have originated in several 
ways. Either, rock of lesser density has been 
mechanically forced into the underlying strata, 
or the underlying strata have lost their heavy 
constituents, or locally mass has been mechani- 
cally removed (for instance by erosion, tectonic 
denudation, fissuring). 

In all these interpretations the physical 
result is the same; in the fairly general 
accepted picture of a floating solid crust of 
several tens of kilometers thickness on a 
substratum behaving like a fluid, this means 
that an upward force is called into existence 
which tries to readjust the disturbed isostatic 
equilibrium. 


Supposed root 


tables for regional isostatic reduction. He 
imagined a sharply limited mass placed on 
the crust; the resulting sagging attains a 
diameter that is larger than the diameter of 
the load, in other words, the compensation 
is more regional than the primary topographic 
element. 

It can also be imagined that a load exerts 
from beneath the crust an upward force on it. 
In this case too, the bulging will have a 
larger diameter than the load, but now the 
topography is the compensation of the 
disturbing mass, in other words, the secondary 
topographic element is more regional than 
the bearing root (see fig. 11). 

Hertz computed that a floating plate 
under the influence of a load P concentrated 
in a point attains a down-bending at this 
point 

P 


aan 
8 Kl? (4), 


Topographie bulging 


+ eu 
a 
ng ae 


Line of compensafron 


Fig. 11 — Schematic section of a crustal bulging caused by an upward force. 


The physical investigations can now be 

divided into three fields: 

l. the investigation into the deformations of 
the crust 

2. discussion on the stresses occurring in the 
erust and their admissibility, and 

3. suppositions about the reaction of the 
substratum. 

Deformations of the: Earth’s crust 

under the influence of a vertical load 


Let us first confine our attention to an 
elastic crust of from 15 to 25 km thickness 
and neglect for the present the influence of 
the substratum, that is, we assume the 
substratum to be a fluid without (great) 
viscosity. 

For this physical configuration, a loaded 
floating elastic plate, Hertz (1884) gave & 
solution (see also Föppl, 1907, pp. 112-130). 
This solution was also used by Vening 
Meinesz (1939) for the composition of his 


in which K denotes the upward pressure 
and l a length defined as follows 


=, 


in which m= Poisson’s constant 
E= Young’s modulus 
and h= the thickness of the crust. 


This down-bending becomes zero at a 
distance of 3.887 X 1, after which it continues 
wavelike, hence with alternating sign, with 
a wavelensth „x[/ y2 and with quickly 
decreasing amplitude. 

A physical”load can now be realized by a 
series of point-loads. The total deformation 
can be obtained by adding the separate 
deformations. 

we simplify our computations — as Vening 
Meinesz did — by neglecting the changing 
of sign and by taking the zero-point of the 


m? E #3 
(m? — 1) 12 K 


(5), 


deformation at a distance of 2.905 Re lstrom 
the centre. This gives the advantage that the 
curve already integrated for the 2-dimensional 
case (see Part I) can be used. The amount 
of the compensating mass, that is, the 
amplitude of the down-bending, was there 
determined by starting from the (justified) 
thought of complete isostatic compensation. 
Here tlie problem will be considered more 
fundamentally and we will start from for- 
mula (4). This formula itself must also be 
integrated for the 2-dimensional case. From 
this integration we get the formula 
cP 

IS grRR 

(6), 
in which c = £(l).>5 

We now have to determine the numerical 
values of Il and f. In Daly (1933) we find 
for Young’s modulus of different rocks under 
normal conditions values varying from about 
1 tot 10 x 1011 dynes/cm? and for m values 
varying from 4 to 10. These results were 
confirmed by Birch and Bancroft (1939). 
Bridgman (1928 and 1929) computed from 
experimental results values for E for arbitrary 
material not of direct importance to geology 
at a pressure of 10.000 kg/cm? (that is the 
pressure which at a density 3 in the earth will 
occur at a depth of 33 km). The changes in E 
were restricted to some percentages of the 
value under normal circumstances. Also from 
experiments Birch and Bancroft comp- 
uted E for various rocks at a pressure of 
4.000 kg/cm? and arrived at values of from 
Be 2173,57, 104 for Sranite-8,4,%x „104, 
Finally Birch and Dow (1936) measured 
the effect of temperature on E for various 
kinds of glasses up to 400° C and stated 
inerements of E of some percentages. Lees, 
Andrews and Shave (1924) established 
somewhat higher values, not greater, however, 
than 10 % more than the normal values. 

If we take for E values of from 1t0 10 x 10! 
and introduce these values into our comput- 
ation, the effect of pressure and temperature 
will most likely be implied. 

With regard to Poisson’s constant, this also 
varies a little as appears from the results 
stated in the literature. The effect of these 
variations of m, however, is hardly appreciable 


5 Vening Meinesz (1946) gives a formula 
for the downbending in the 2-dimensional case 
from which it appears that ce = £f(l) = 2.21, This 
formula gives values that only deviate about 
1% from the values that were obtained graphi- 
cally. 
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in (5), so that for m 


4 
in the factor % m? 
m?—| 
the value 4 may be taken. 

K denotes the upward force that corres- 
ponds to the density of the substratum minus 
the density of the seawater, as the upheaval 
in question is still beneath sea-level. There- 
fore we have taken 2 for the value of K. 

For h we shall take 15, 20 and 25 km 
successively. 


T'hus I varies from 22 to 55’km. For h = 
15 km, Il becomes 22-33 km, for h = 20 km, 
l becomes about 44 km, and for h = 25 km, 
l is about 55 km. We will adopt for | values 
of 22, 33, 44 and 55 km. 

The 2-dimensinal disturbing mass has been 
taken firstly 60 km wide and 20 km deep, 
secondly 120 km wide and 10, 15 and 20 km 
deep successively, in all cases with a density 
difference 0.6. From this f follows (formula 6). 
Combined with 1 we thus obtained 16 bending 
curves. The shaded part in fig. 11 represents 
the topographic bulging, i.e. the topographic 
compensation of the root. It is fairly wide 
(from 150 to 400 km). This bulging could 
perhaps be identified with the broad ridge 
that appears on the oceanographic maps south 
of Java at the place of the Vening Meinesz 
zone (see also profile VI and VII). Naturally 
the bulging cannot be separated exactly from 
the rest of the topography, but at any rate, 
it is not wider than about 250 km. Applying 
regional isostatic reduction, the ridge is 
supposed to be compensated according to the 
broken line in fig. 11. The distribution of mass 
supposed in regional isostatic compensation, 
and as follows from the picture assumed here 
— a bulging caused by a root — now differ. 
A negative mass in the centre and positive 
masses on either side are causing gravity 
anomalies. The resulting anomaly picture was 
computed. 

Tabel IV presents the results of these com- 
putations. Since the assumption of a root of 
120 km by 15 km takes but an intermediate 
place amongst the others, the pertaining 
results are left out. 

Briefly summing up, the following may be 
said. The computed anomaly picture, pertaining 
to the assumption of root dimensions 60 by 
20 km was too narrow with respect to the 
negative anomalies, the anomalies resulting 
from a root of 120 km by 10 and 120 by 15 km 
were too small. 

In case of 1 = 22 km and root dimensions 
of 120 by 20 km, the upheaval appeared far 
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TABLE IV — ISOSTATIC ANOMALIES (MGALS) RESULTING FROM CONSTRUCTIONS (SEE TEXT). 


u ET u u u EZ EEE DEUTET 
Distance from centre in km 


Regiona- Value Root dimensions 
lity of in km 0 30 60 90 120 150 180 
= 0 20 —110 — 61 6 31 26 18 11 
A 210 E69 — 56 —17 21 28 18 13 
120 x 20 —126 —111 —39 38 49 39 29 
44 60 x 20 —118 — 41 20 36 27 11 4 
120 x 10 — 66 55 —11 21 30 19 10 
208x220 —t31 —111 —92 50 52 38 23 
33 60 x 20 150 — 55 32 37 25 15 5 
120 x 10 —ı68 —ı59 —6 31 26 15 8 
120. x 20 —127 —110 —17 65 54 36 18 
22 60 x 20 —102 — 40 42 50 31 18 7 
202x210 — 60 — 54 —3 31 26 15 gi 
1209820 —-115 —105 —12 60 53 30 14 
R=0km 44 1207220 — 83 — 71 —5 46 36 17 4 
33 120 x 20 — 64 — 60 3 47 21 5 2 

rounded shape 

44 of root — 51 — 48 —28 4 32 26 13 


too high as compared with the ridge south of 
Java (more than 5 km in the centre), in case 
of 1 = 55 km and the same root dimensions 
this upheaval is far too broad (more than 
350 km). 

A study of fig. 11 tells us that the anomalies 
have to become smaller if we apply a local 
isostatic reduction for the ridge. To some 
extent this happens to be in accordance with 
the facts. In table V are given the isostatic 
anomalies for Vening Meinesz’s stations 
179 and 183 as resulting from the various 
methods applied for the isostatic reduction, 
followed by the values we calculated for the 
gravimetric effect above the centre of the 
upheaval of the difference between a local 
compensated ridge and a ridge caused by a 
root of 120 by 20 km if 1 = 33 and 44 km 
(see table IV). Both stations are situated at 
the extreme minimum of the anomaly curve 
and so this comparison seems to be justified. 

From the table it is clear that the decrease 


of the anomalies is far too great if 1 = 33 km. 
In the other case the construction may still 
represent an approximation of. the actual 
situation. 

Let us now compare the regional isostatic 
anomaly curve of the remaining construction 
with the gravity profiles VI and VII for the 
part south of Java. For this purpose the 
computed anomaly curve is indicated schemati- 
cally in Profile VI and VII. It immediately 
strikes the attention that the positive values 
resulting from the construction too closely 
approach to the negative anomalies. This 
discrepancy cannot be forestalled. If the 
thickness of the crust were taken twice as 
large for instance, then I would become 1.6 
times larger, that is to say, l becomes 88 km 
instead of 55 km. From the values in table IV 
it is clearly seen that the outward shifting of 
the positive anomalies do not keep pace with 
incerement of I. Furthermore, the gradient of 
the negative part of the curve becomes too 


TABLE V — RELATION BETWEEN THE DEGREE OF REGIONALITY OF THE ISOSTATIC REDUCTION 
AND THE ANOMALY AT THE CENTRE OF DISTURBANCE. 


Regionality in km 


—127 


| —134 


Vening Meinesz’s 
Station 179 u 
1dem 183 T='30km —126 
T = 20km —115 
I = 44 km eg 
naar Er 
Idem en 
es in =.lökm | 84 


small. In fact, in this wäy a sufficient distance 
between the positive anomalies cannot be 
reached. Moreover, if I is larger, the diameter 
of the bulging is no longer in agreement with 
the actual state of things. It is therefore fairly 
certain that in this manner no solution of the 
gravity field can be found. 

The value of these conclusions might still 
be doubted by objecting that a rectangular 
shape of the root — as assumed — physically 
- speaking is quite improbable. For this reason, 
yet another construction was computed with 
Il = 44 km and a more rounded shape of the 
root. In this case the negative anomaly proved 
to diminish considerably (up to —51 mgals 
in the centre, see table IV), although the total 
ot the disturbing mass had already been taken 
somewhat larger than in the case of a 
rectangular root of 120 by 20 km. Should the 
total of the disturbing mass be enlarged even 
-more, then the central rising would become 
improbably large. 


The strength of the Earth’s crust 


We shall now try to examine what stresses 
occeur in the earth’s crust when such a load is 
applied and whether the crust is able to sustain 
these stresses. Hereby we shall continue 
assuming that the effect of the substratum is 
negligible. 

Theoretically the problem is not so simple. 
In Hertz’s solution the stress distribution in 
the plate is not given. Moreover, there is no 
sufficiently adequate theory on the conditions 
causing failure. The existing maximum stress 
and strain theories are all extrapolations of 
results secured from yielding tests under 
restrieted stress conditions. Bridgman (1912) 
succeeded in giving a clear experimental 
illustration of this. 

Even if the problem is approached from the 
angle of the lattice energy one is up against 
difficulties not yet overcome. 

The best thing to be done in this situation 
is to see whether experiments have been 
carried out concerning the conditions of failure 
under analogous circumstances. 


Bridgman’s work (1935, 1936, 1937) is an 
excellent guide here. He loaded material in 
an apparatus, the schematic section of which 
is reproduced in fig. 12. Bridgman calculated 
that the pressure between the pistons B and 
the anvil A approximately becomes of hydro- 
static character as a result of the deformation 
of the pistons at the enormous pressures he 
applied. By measuring the force needed to 
turn the handle C, information was obtained 
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about the shearing strength. The applied 
pressure mounted up to 50.000 kg/cm?, i.e. the 
pressure that occurs in the earth’s crust at a 
depth of 150 km if the density is 3.3. 

Among the many examined substances there 
were four minerals. The results of the experi- 
ments on these minerals are of direct im- 
portance to geology. They are 


mullite (synthetic) — 3 Als032SiOs 
sillimanite — Als03Si03 
andalusite — AlsO3SiOs 


and lepidocrosite — Fe&0;3.Hs0, 


Fig. 12° — 
apparatus. 


Schematic section of Bridgman’s 


Of these andalusite is most representative of 
the rockforming silicates. Figure 13 shows 
Bridgman’s results. The value of the shearing 
strength at a pressure of 1 kg/cm? has 
been derived from Daly (1933) and is the 
shearing strength of granite. For comparison 
also the curves of aluminium and of paraffine 
are given, both representing extremes: alumi- 
nium showing a fairly straight curve with a 
factorial increase of 10 at 50.000 kg/cm?, 
paraffine showing a steeply rising curve with 
a factorial increase of 900 at 30.000 kg/cm?. 

At the moment we are only interested in 
the part of the curves up to 10.000 kg/cm?, 
that is, to a corresponding depth of 30 km. 
For this purpose we have connected the value 
of the shearing strength of granite at atmos- 
pherical conditions with the mean of the 
strength of the four minerals at 10.000 kg/em?. 
This part of the curve can be approached with 
the formula 

un = 0.17 p + 300) kg/cm? 
(7) 


in which 
7. >= the shearing strength and 
p the hydrostatic pressure. 


Let us now consider the shearing stresses 
that will occur along the plane AB Intel® 
In the case of a root of 120 by 20 km at a 
value of 1 = 44 km and with a thickness of 
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the crust of 20 km — the construction that 
gave the best gravimetric results with regard 
to the negative anomalies — in this plane a 
shearing stress of 1.200 kg/cm? has to be 
sustained. And if we take the thickness of the 


15000 
kg/cm! 


10000 


SHEARINO STRENGTH 


5000 


2000 


on the shearing strength has been examined; 
an examination, moreover, only in connection 
with a plane in which indeed very large 
shearing stresses are to be expected, but in 
which probably not the largest shearing 


= 
20000 


= 
50000 
kg/emt 


u Be 
30000 40000 


HYDROSTATIC PRESSURE 


Fig. 13 — Diagram showing the relation between hydrostatic pressure and the shearing-strength 


of several substances. 


crust at 25 km, the figure becomes 960 kg/cm?. 
Both cases are represented in fig. 13, with two 
horizontal lines. Supposing that there is 2.5 km 
of seawater upon the crust and assuming that 
the crust’s density is 2.7, the hydrostatic 
pressure ranges in the first case 

of from 250 kg/cm? to 5650 kg/cm? and the 
shearing strength 

of from 340 kg/cm? to 1260 kg/cm?, 
in the second case the hydrostatic 
ranges 


of from ‘250 kg/cm? to 6900 kg/cm? and the 
shearing strength 


of from 340 kg/cm? to 1490 kg/cm?, 


Thus, in the first case the strength is 
evidently far too small; it is rather small in 
the second case, and probably still insufficient. 


So far only the effect of hydrostatie pressure 


pressure 


stresses do arise, 
stress distribution. 


The effect of time, (higher) temperature 
and inhomogeneities on the strength of the 
material must not be neglected, however, This 
influence is a weakening one. The effect of a 
higher temperature is obvious, the weakening 
influence of inhomogeneities on the material 
is also known. As to the effect of time, 
Griggs (1936) performed interesting experi- 
ments, He arrived at the conclusion that the 
elastic limit was not influenced by the 
duration of loading, but that the ultimate 
strength was. The ultimate strength approaches 
asymptotically to a time-independent value. 

In this matter it is difficult and dangerous 
to make quantitative pronouncements. We 
may predict, however, that the combination 
of all factors render it very unlikely that a 


considering the intricate 


crust of 25 km thickness has a sufficient 
‚strength and it may be asserted that a crust 
of 20 km thickness is not strong enough. 


Consequently, if no other forces are assumed, 
the rate of upheaval of the Vening Meinesz 
zone would be determined by the rate of 
reaction of the substratum to disturbances of 
the isostatic equilibrium. 


The reaction of the substratum 


Up till now we have entirely neglected the 
influence of the substratum. We shall now 
investigate what part is played by it. 

First of all we may suppose that the 

substratum is sufficiently strong to sustain 
elastically the stresses sketched in the fore- 
going. We ihen could apply any anomalous 
distribution of mass whilst no reaction of the 
substratum worth mentioning sets in and in 
this way any solution of the gravity field 
might be put forward. 
“ However, the phenomenon of isostasy 
indicates a rather small strength of the sub- 
 stratum and even if the adjustment of isostatic 
equilibrium has to be described to very slow 
processes like creep, there is still the con- 
vincing instance of the postglacial upheaval 
of formerly ice-capped areas, Fennoscandia 
and North America. 


Jeffreys (1952, p. 326 ff.) disputed the 
value of the Scandinavian phenomena as an 
argument for isostasy basing this mainly on 
the absence of uplift of the region of the 
Hudson Bay in Canada. Gutenberg (1941), 
however, demonstrated clearly enough that 
an uplift of this region is actually taking 
place. But even when a recent uplift of this 
area could not be shown, it need not surprise 
us. For it is fairly certain that the ice cover 
of such a vast area as the North American 
continent had a more complicate form than 
the Fennoscandian ice cover. Considering 
fig. 1 of Vening Meinesz’s latest paper on 
this subject (1954), it is obvious that besides 
uplift also subsidence occurs; in case of a 
somewhat more complicated form of the ice- 
sheet, it may happen that a more limited 
area remains at the same level for a certain 
time. At all events it is beyond dispute that 
the Canadian Lakes are being tilted, while 
another explanation is hard to find that sticks 
so closely to the facts as adjustment of 
isostatic equilibrium. 

A detailed gravimetric survey of the North 
American region has not yet been carried 
out; the only observations that have been 
published, however, do indicate the negative 
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anomaly of the area to be expected there 
(Innes, 1948), 


In the NW European area of glaciation a 
striking agreement between theory and obser- 
vations can be established (see e.g. Vening 
Meinesz, 1954). In the following pages we 
shall confine ourselves to the more detailed 
data on the upheaval of Fennoscandia and 
extrapolate the results to Indonesia. So we 
assume that the zone of weakness, the 
„asthenosphere’”, that appears from the post- 
glacial upheaval beneath NW Europe and 
North America, is present all round the earth. 


For a mathematical treatment of the problem 
two ways are open to us, each of them based 
on different physical assumptions, namely, 
plastic and viscous flow. 


In the publications known to us a solution 
of the problem was always founded on the 
conception of viscous flow, that is, the sub- 
stratum is without any strength, although 
most authors are aware that this implies an 
approximation. So whenever the term visco- 
sity is met with in the following pages, we 
had better read pseudo-viscosity. 


Of the calculations carried out on the 
postglacial upheaval of Fennoscandia, perhaps 
the most important are those of Van Bemme- 
len and Berlage (1934), Haskell (1935), Vening 
Meinesz (1937 and 1954), Niskanen (1939 and 
1948) and Gutenberg (1941). 


Van Bemmelen and Berlage 6 concluded to a 
kinematic viscosity coefficient 0.3 x 1020 with 
a thickness of the viscous shell of 100 km, 
Haskell got 3 x 1021, Vening Meinesz (1937) 
obtained 9 x 1021 and (1954) 2.5 x 1021, 
Niskanen arrived at 1.1 x 1022 and Gutenberg 
got 7 x 1021 c.g.s. units. 

Niskanen’s more recent publication is an 
elaboration of the problem with the aid of 
subsidiary assumptions. 

The result of Van Bemmelen and Berlage 
is surprisingly small compared with those of 
the others. To begin with the latter will be 
regarded as a group, next the difference 
between these results and Van Bemmelen and 
Berlage’s will be dealt with. 

By comprising the various formulae into the 
form 

N °} 
v”T2kb 
(8), 


6 Their paper stated a normal viscosity coef- 
ficient of 13 x 1020, This should read 0.8 x 1020, 
which figure divided by the density 3 gives the 
kinematic viscosity coefficient. 
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TABLE VI — COMPARISON OF THE VARIOUS METHODS OF DETERMINATION OF THE K 
MATIC VISCOSITY COEFFICIENT WITH THE AID OF FORMULA (8). 


Name of author Meaning of k Meaning of b Value of k Value of b 


quantity inversely propor- 


a a ee. 


Vr 


= ra ä _12 BB 
Haskell "vertical dimension of | the disturbance and depen- |+ Yn23x 10 /sec|+ 1.1 x 10-75 
disturbance | dant on *k - 
Vening SUN DEIE r/diameter of the distur- 


vertical dimension of dis- 23 x 10-'?/sec 


turbance 


22 x 10-%a 


Meinesz (1937) bance 


damping constant depen- 
dent on the time elapsed 


: kell’ 1? /sec 13 x 10—°/e 
seien since the beginning of the en 3.4 x 10 / 
upheaval u. 
Gutenberg idem idem 7x 10-!?/sec| 13 x 10-3/e 
Dee ash idem idem 71x 10-1?/see| 33 x 10-%/& 
einesz 
Van Bemmelen a DL : = 3L? Formula is only indirectly used fo 
d Berlage vertical of dimension dis-| b = —— R determination of the viscosity 
er turbance 2 h’n? 
22 x 10-°/en 
id 1/f (L/zu, h), see text SD — 12 a —— ee 
Extended formula —ı F(L/n. h) 23 x 10 e_ 29 x 10-Ya 
m 
in which meters for elastice rebound is taken into 
» = the kinematic viscosity coef- dccount. Gutenberg and WNiskanen found 
fieient ‘ their value for b in a somewhat different 
= the value of gravity, while way from the data than Vening Meinesz; 
gI 7 ’ 


the difference between the outcome of 
Gutenberg and Vening Meinesz is resulting 
from this. Since it is not possible to prefer 
the one method of determination of b to the 
other, we cannot go further than establishing 
that a value for „» between 2 and 7 x 1021 
is probable. 


k and b are constants, a good synopsis can be 
obtained of the theoretical approximations 
made use of and of the adaptations to the 
data. The various values of b and k according 
to the various solutions are represented in 
tabular form (Table VI). 


The theoretically most elaborated solution 


is the one by Niskanen. In determining the 
viscosity coefficient he started from a total 
amount of upheaval of 450 meters. The data 
before the year 6800 B.C. he found too vague 
to be introduced into his calculations. That 
total amount, however, is of great importance 
for the determination of the damping constant 
k (see Vening Meinesz, 1954, formula (8)). 
Thus Niskanen arrived at too great a 
value for y», Gutenberg (1941) and Vening 
Meinesz (1954) followed both the same way 
to determine the part of the disturbance 
dispersed up to 6800 B.C. The latter argued 
that this amount must have been 270 meters, 
so that the total amount of the disturbance 
in 7700 B.C. must have been 620 meters, if 100 


Haskell’s value does not differ much from 
it, but Vening Meinesz (1937) demonstrated 
that on account of the neglecting or underesti- 
mating of the amount of upheaval to be 
expected, an error crept in, Hence, Haskell’s 
result has to be 9 x 1021, Instead of Haskell’s 
the corrected values for vo, k and b, are 
presented in table VI. 


That Haskell’s result does after all differ 
a little from that of Vening Meinesz and 
Gutenberg, must be regarded as a result 
ot the velocity distribution, assumed by 
Haskell to simplify his solution. Owing to 
this his formula is but a first approximation. 
Another cause of the difference in result is the 
different value for b (see table VI), a conse- 
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TABLE VI — Continued. 


Depth of ? ee : : f 
E..0 y in c.g.s. units of difference with Vening Meinesz |Half-life for the Vening 
J Meinesz zone 
infinite 9x 101 simplifying velocity distribution from 
which the formula is but an approximation 
infinite 9x 10? simplifying velocity distribution 55.000 years 
7%) different adaptation to data from which 
infinite 11x 102 the total time of upheaval is different; 
() other way of determination from the 
data (see text). 
Infinite - 7x 1021 (6) other way of determination from the 
data 
| infinite 2.5 x 102! < 10.000 years 
100 km Ile other assumptions about the depth of the 
20 km 2x 107 substratum ee 
infinite 9 x 10°! simplifying velocity distribution 55,000 years 
20 km 7x 10" slight depth of the substratum 275 years 
quence of the circumstance that Haskell Yu ou Ju du 
considered the upheaval along his longer axis R u d% is %y Te == 
(over Oslo). The difference between Vening 
Meinesz’s results (1937 and 1954) is a (> R o?u ) ar 
consequence of the simplification made in the ,} I.x 9Yy° 92? oe 9x 
former paper, in which the disturbance was | j 
represented in the form of a single sinus gu u 2 eu ge ee 
function. The forces occuring now are as a of op 9Yy 92 
matter of fact greater than in the case of eo ey 32y 13 
single disturbancee and consequently the _» ( +: ; = I) = vera) 
resulting viscosity coefficient is also greater. IX 9Y 92 0 0Y 
So there remains the factorial difference of 9w ar aw 1. yw I 3w 
70 between the results of Van Bemmelen 9 Ix ey 92 
and Berlage on the one hand and the other 2 2 1 
authors on the other. We will now try to iind: =& 5 ( Pr 0 a: g 2) = Ko Ip 
the cause of this difference. If it is assumed, dx? aoy” ie 0 09z 
like all the above authors did, that the fluid (10) 


is incompressible, that is, 


(9), 


the equations of motion for viscous fiow 
become: 


the velocity compo- 
ponents 

the kinematical 
viscosity 

the volume forces 


— the density 
— the pressure (Web- 


1949, p. 550). 


ster, 
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By assuming the problem to be 2-dimen- 
sional, as all did, except Haskell, and in- 
serting the gravity as the only volume force, 


the equations become 


u du ou EN 
nie) 
1lo9p 
re: 
a (= a) 
ot 9x 22 Pr 
l op 
0 9 
ny* 


In order to come to a solution, further 
must be made. 
based their 
equations on Defant (1932), who derived a 


simplifying assumptions 


Van Bemmelen and Berlage 


solution especially for the problem of sea- 
waves, He neglected the height of the waves 
with respect to their length and thus supposed 
horizontal laminar flow, that is, all terms in 
w become zero. From this it follows since the 
fluid is considered to be incompressible (for- 


au 9°u 
mula 9) that —— and <=, become zero, In 
0129 9x 


ß) 
addition he assumed that the term ss may 
[e 


be substituted by the horizontal pressure 
gradient of the disturbing mass 


RI@ 
eg SE 
By this the equations reduce to 
au Urne 
ot 9 x 09x” \ ä 
R (12) 


and in order to come to a solution Defant 


gave the equation of eontinuity 
ordu dc 
ee zZ olS) —— 0) 
Zn 0 ot 
(13), 


in which R = the thickness of the pseudo- 
viscous layer. 


Van Bemmelen and Berlage obtained from 
this the formula for the rising in the centre 


—1,g hin? 0, 
al 


a - a 


(14), 


in which L represents the diameter of the 
disturbance. (This formula is used only in- 
directly for the determination of the viscosity, 
but this does not affect the hydrodynamic 
principles.) 

The assumption of laminar flow involves 
that the flow may occur only down to a 
rather slight depth. If this depth is taken 
too large a vertical flow must occur in the 
This 
vertical flow is no longer to be neglected. 


upper layers to feed the lower ones. 


So the results obtained by Van Bemmelen 
and Berlage for the viscosity coefficient by 
means of their formula on page 43 (op.cit.) 
cannot be compared with the results of the 
other authors, who all assumed the layer to be 
infinitely thick. For when Rh in this formula 
is enlarged, in other words when the layer 
is made thicker the boundary conditions are 
disregarded. 

Vening Meinesz (1937) and Haskell supposed 
that the time derivatives and the first space 
derivatives can be neglected, that is, they 
assumed very small velocities, quite a different 
simplification than the laminar flow of Defant. 

Hence we get the equations of motion of 
Vening Meinesz 


oa 92 09 
ee EEE 
9x2 l 2° 0 zZ g= 


(15). 
(It has already been said that Haskell started 
from the 3-dimensional case, So he had the 
corresponding equations of motion.) 


By assuming a periodie disturbance Vening 
Meinesz obtained the formula 


(16), 


In this formula Ah, the thickness of the layer, 
is no longer found, the phenomena take place 


Liu 


ARE 1 u 


ne u nn nun u re en u 


7 


N ER % 
> 


in an infinite thick layer. 7 

When the formula of Van Bemmelen and 
Berlage is compared to the one of Vening 
Meinesz, then it strikes the attention that L, 
the diameter of the disturbance, occurs in 
the first formula (14) to the power —2 and 
in the second (16) to the power +1 in the 
e-exponent, and further as has been said 
and explained already, that enlargement of h 
in Van Bemmelen and Berlage’s formula gives 
no real results. £ 

Now it is very important to know how 
these formulae have to be extrapolated, 
formula (14) to the depth, formula (16) in 
assuming finite depths. An extension of the 
formula of Vening Meinesz will therefore 
be given below, an extension that enables a 
further analysis. 

For this purpose we start from the equations 
given by Vening Meinesz (1937). Equations 


- (15) here formed the equations of motion, the 


condition of incompressibility was expressed 
by formula (9), the stress equations were 
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that u and w 
in 


Vening Meinesz assumed 


are functions of z, multiplied by e i or in 


other words that u and w are periodic in the 
direction of the X-axis with a half wave 
length L. Through this the equations obtain 
an other form. By elimination follows the 
equation in w 

nn? 9°w zı* 3 
—'2 3 rer 


EN 
el 
(18), 


in which we may leave out the last term, since 
the general component of the gravity does 
not actually contribute to the flow. The 
general solution is then 


7 -77 
WEoNeLııceL. 
TU 77 
er. 1 en 
(19). 
Vening Meinesz now assumed as boun- 


9w SR 
u 2vo = dary conditions that u and w become zero 
n = for z = infinite and that 'y = 0 for z = 0 
9 gw We, on the other hand, will assume that u 
AI Se and w are zero for z = h and further that 
%y — 0 ifz — 0. So the constants A, B, C and 
au aw Dee 
z =mne|l + — D can be determined. 
y 920% The ratio L/r in the coefficient of (16) 
(17), appears to be substituted by a very compli- 
cated function of [/n andh. 
in which Gx and 0, = the normal stress 
components So 
Dh 
z, > the shearing stress eh IE 
component en 2n t 5 
and p = the pressure = —%(o. + 0;z) (20) 
in which 
BR N 
EN ER IE L EURE 
r(E " it a a 
n Du TU TU g7 en 
> h E h h 
3 75 ERDE, L Hi L 
we 2=e se ne 
IT IT“ 3 
rl (20)a. 


7 Niskanen made the assumption that only 
the first space derivatives of the velocity are 
negligible. Moreover he did not start from a 
simple sinusoidal disturbance of the equilibrium 
but he approached the rising curve with Fourier’s 


series. Theoretically this solution is more elegant. 
As appeared in the foregoing (table VD, the 
results do not greatly differ and furthermore the 
differences could for the greater part be traced 
back to other auxiliarly assumptions in the 
treatment of the problem. 
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By numerical calculation it then appears 
that in case of a small h, L becomes very 
important in the denominator of the e-expo- 
nent and occurs nearly as a square. The 
formula approaches in this region the so- 
lution of Van Bemmelen and Berlage. For 
infinite h the function reduces to |/n and 
the solution becomes identical- with formula 
(16). The function being continuous, in case 
of a gradually increasing h, L has to become 
more important in the enumerator, 

So this formula fills the gap between (14) 
and (16). 

The question may be asked what the change 
of the exponent of L physically means. The 
derivative of the general formula (20) repre- 
sents the velocity of decrease of the distur- 


bance, 
dd BORN Le 
ae 2 rl. IE 


(21). 


We now read:.at a greater depth the velocity 
of decrease is proportional to the wave length, 
at smaller depths this velocity is inversely 
proportional to the square of the wave length 
and further dependent on h. In the latter case 
the friction with the bottom of the vessel, i.e. 
the upper surface of the solid layer beneath 
the viscous one seems to act considerably 
retarding. 

In the hydrodynamics of non-viscous fluids 
we find an analogy of this. The following 
phenomenon occurs in migrating waves: at 
greater depth of the basin the velocity is 
proportional to the wave length, at smaller 
depth the velocity is independent of the wave 
length and proportional to the depth. So with 
regard to the smaller depth there is still some 
difference. However, it can easily be imagined 
that the change of exponent of L in viscous 
fluids is even greater. 

We shall now try to extrapolate to Indonesia 
the results to be obtained by means of the 
preceding formulae. 

In assuming an infinite thick layer of sub- 
stratum Vening Meinesz (1937) gave for 
the kinematie viscosity coefficient 9 x 1021 
c.g.s. units. 

Assuming a thin viscous layer we will have 
to start from the extended formula; at a thick- 
ness of the layer of 20 km, » becomes 7 x 1017, 
If, with the same assumption we compute y 
with the aid of Van Bemmelen and Ber- 
lage’s formula, then y turns out to be 
2 x 1017, which is a fairly good agreement. 


The remaining difference will be a conse- 
quence of the fact that Van Bemmelen and 
Berlage inserted a value for the diameter of 
the disturbance of 1200 km, whereas Vening 
Meinesz took 1400 km, and further of the 
circumstance that Van Bemmelen and Berlage 
based their computations only indirectly on 
formula (14). For, inserting the e-exponent 
K, computed by Vening Meinesz in formula 
(14), the viscosity increases by a factor 3 and 
thus becomes about equal to the value ac- 
cording to formula (20). 

Supposing that these two assumptions — 
an infinite thick layer and a layer of 20 km 
thickness — represent extremes, the conclusion 
may be made that the value of the kinematic 
(pseudo)viscosity coefficient lies between 
9 x 1021 and 2 x 1017, 

By means of formula (14), (16) and (20) 
it will now be tried to compute a half-life for 
the anomalies of the Vening Meinesz zone, 
that is, the time in which the anomalies will 
be halved. The width of the Vening Meinesz 
zone is taken 200 km for our purpose. 


The general formula is 


nei 
seu® 
(22), 
from which follows the half-life T 
In2 0.69 315 
sr K 
(23) 


For an infinitely thick layer K — eg en 
? E14 
which gives with el 
o 
0.4 x 10-12, from which follows T = 55.000 
years. 


=: 9 x 1021 a value 


For a finitely thick layer K= a f ee h). 
205: 


Fork = 20 kmand » = 7x 1017 this gives 
8.0 x 10-11 and T-becomes 275 years. 
In the formula of Van Bemmelen and 


Is gh’a?o 
9.6” 


v= 2x 1017 for T a value of 66 years. 
Considering further the rate of rising in 
the light of Vening Meinesz’s most recent 
publication (1954), then from fig. 1 (op. eit.) 
it clearly follows that the half-life for Indo- 
nesia will be smaller than that of Fennoscandia, 
For a relatively narrow and high disturbance 


Berlage K = which gives with 


spreads out more quickly than a wide and 
less high disturbance. 

The various half-lives are also represented 
in table VI. 

Another way to approach the upheaval of 
Fennoscandia is by assuming that the sub- 
stratum indeed has a definite strength and 
that the flow is of plastic character. Unfortuna- 
tely there are no equations to solve this problem 
in an exact way. Nevertheless attempts will 
be made to arrive at some conclusions about 
the substratum by means of qualitative 
considerations. There are fundamental dif- 
ferences in treating the substratum as a 
viscous liquid and as a plastically behaving 
solid, as may be shown in the next few lines. 

In the case of viscous flow the stresses are 
propertional to the displacement velocities. 
In the case of plastic flow, however, only a 
certain amount of stress can be taken up, 
determined by the Mises-Hencky flow 
- condition (see Reiner, 1949, p. 111). The 
surplus of mass-force causes an acceleration, 
so flow occurs with increasing velocity. The 
occurring stresses are therefore independent 
of the displacements. 3 

To illustrate this contrast, fig. 14 and 15 
have been taken over from Reiner’s book. 
Figure 14 represents viscous flow, figure 15 


PERFORATED 
PISTON 


LIQUID 


Fig. 14 — Model for viscous flow (after Reiner, 
1949). 


FRICTION II 


Fig. 15 — Model for plastic flow (after Reiner, 
1949). 


plastic flow. The latter clearly demonstrates 
the analogy between plastic flow and the 
movement of a body which has overcome a 
certain friction. 

This representation of plastic flow is greatly 
idealized and occurs under normal conditions 


8 See also. J. M. Burgers “Mechanical consider- 
ations ete.” in the First Report on Viscosity and 
Plastieity (Proc. Kon. Acad. Wetensch. A’dam, 
1st Section, Vol. XV no. 3, 1935). 


only just after the yielding-point of the 
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material. For the deformations setting in after 
this, again mathematic formulations can be 
put forward, which involve a combination of 
elastic, plastic and viscous deformation. 

From Bridgman’s experiments (1951) it 
appeared, however, that this idealized formu- 
lation of plastic flow can very well be used 
in the case of plastic flow under high pressure: 
”But at high pressures and very large defor- 
mations the ideal state is found to-prevail over 
an indefinitely wide range of information”. 
So it is fair to start from the ideal plasticity 
and to put aside the other assumptions. 


Imagining this flow to occur in an infinitely 
thick layer, then the Mises-Hencky flow 
condition gives us the boundary conditions of 
our problem; at a certain distance from the 
load given shape by the plane „, another 
plane (7) is formed, for which u and w, the 
velocity components are zero. This plane is 
a function of the load and of the Mises-Hencky 
flow condition. 


y= (9,950) 

(24), 
in which os expresses this flow condition. 
Because of the. decrease of the disturbance the 
plane will move towards the centre of the 
disturbance and thus reduce the plastic region. 
This reduction may be so important that the 
accelerating influence is completely abolished 
and the movement may even be retarded. 
This would be true for Fennoscandia, where, 
as has been said, the upheaval can be repre- 
sented by the expression 


w—-Kc ‚from which BEL SRT 
(22), 
K is now a function of yp , So that 
< ud 2 $ ar ), ‚ 
Knete Fmoslt 
(@BBIE 


With this formula no exact computations 
can be carried out yet. In the following we 
have to confine ourselves to qualitative con- 
siderations. 

The plastic flow can be dissolved into two 
components, a spatialand a local one (fig. 16a). 
If these two are separated, the effect of the 
spatial component may be represented as in 
fig. 16b, the effect of the local component as 
in fig. 16c. It may now be asked what happens 
if the disturbance is made smaller and higher 
without changing the total mass. The spatial 
component will remain the same or diminish 
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a little. whereas the local component will 
remain equal or — what is more probable on 
account of the greater pressure gradient — 
become larger. 

It is now probable that the velocity will 
become greater, but it is at any rate fairly 
certain that it will not become smaller. 


Disfurbence 


N 


TE 
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Fig. 16a — Velocity distribution in case of plastic 
flow of the substratum; A designs the local 
component, B the spatial component. 

b — The effect of the spatial component. 

ce — The effect of the local component. 


The deficieney of mass of the Vening 
Meinesz zone is certainly not smaller than 
the original shortage of mass after the melting 
of the ice cover in Fennoscandia. 

With some caution the expectation may 
now be put forward that when the reaction 
of the substratum has to be regarded as being 
of plastic character, the half-life of the 
anomalies of the outer arc in Indonesia will 
not be greater than the half-life of the distur- 
bance in Fennoscandia, and at the most will 
not exceed 50.000 years either. ® 


Preliminary conclusions 


From these results it may be concluded: 


1. Assuming a thin viscous layer beneath 
the crust sensu stricto, no retardation of 
geological importance is to be expected in the 
adjustment of the isostatie equilibrium in case 
of narrow and elongate disturbances. 

2. When assuming an infinitely thick viscous 


layer, the half-life of such a disturbance is 
at most 55.000 years. 


3. Supposing that the elastic properties of 
the substratum will oppose a viscous behaviour 
so that the problem has to be dealt with 
according to the laws of plastic flow, the half- 
life of a narrow and elongate disturbance 
cannot be more than some tens of thousands 
of years either. 


4. Once isostatic equilibrium is attained, 
remaining irregularities in the mass distri- 
bution cannot account for gravity anomalies 
as established in West Indonesia, as follows 
from the elastic deformation laws. 


5. In addition, a thin crust of 20 to 25 km 
thickness would not be strong enough to 
sustain the great stresses occurring in such a 
configuration. 


FINAL CONSIDERATIONS AND 
CONCLUSIONS 


In the first part it was seen that the intro- 
duction of a geologie and a geologice isostatic 
reduction for the sedimentary basins did not 
entail large changes for the gravity anomalies. 
Important deviations from the picture positive- 
negative-positive could be traced back for 


9 The Subsidence of the North Sea Region. 

Umbgrove (1950) indicated four (possible) 
causes for the present subsidence of the Low 
Countries: A. Rise of the sea-level, B. Compaction 
of sediments, C. Subsidence as reaction to the 
postglacial upheaval of Fennoscandia after the 
melting of the ice and D. a general subsidence of 
the sea floor. 


As to the third cause (C), Van Bemmelen 
and DBerlage (1934) computed this effect and 
arrived at the conclusion that this subsidence is 
no longer very important. Niskanen (1939) 
indicated this subsidence, but did not give any 
figures. Vening Meinesz (1954) carried out some 
computations on this point and concluded that 
the maximal remaining subsidence in ceonsequence 
of the uplift of Fennoscandia will be but a few 
meters. A new aspect in the last-mentioned 
publication is the possibility that the melting of 
the root of the Alps plays an important part in 
the epeirogenetic movements. On this point we 
may refer to the paper in question. In the above 
it was argued that the assumption of a strength of 
the substratum will probably have but little in- 
fluence on the theoretical half-life of the anoma- 
lies of the Vening Meinesz zone. The question 
is different for the boundary regions of the 
disturbance. Ata certain distance from the centre 
the stress differences might become so small that 
plastic flow becomes impossible. 


In connection with the projects for the endiking 
of the sea-arms and the other endikements it 
thus could be of hish importance to carry on the 
investigations into this question and to start then 
from the supplementary assumption of a plastic 
behaviour of the substratum. 
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the greater part to the” geology. As for the 
great uncertainties, no computations were 
carried out on this point. 

In the second part the conclusion was 
arrived at with a fairly high degree of cer- 
tainty that the anomalies of the Vening 
Meinesz zone with the supposed thickness 
of the crust, i.e. 25 km, will have a life-time 
of some tens of thousands of years if no other 
forces are presumed, This time is — geologi- 
cally speaking — far too short, and so the 
conclusion cannot be escaped that this solution 
is impossible., 

This conclusion was reached by assuming 
that the crust behaves as an elastic plate 
floating in a fluid. Once equilibrium is obtained 
after application of a reasonable load, this 
configuration does not give sufficient agree- 
ment with the actual condition with regard 
to its gravimetric effect. In addition, this 
equilibrium is not stable, since in all proba- 
bility the strength of the crust’s material is 
insufficient to sustain the enormous stresses 
occurring in it. The half-life for the adjust- 
ment of the equilibrium was extrapolated from 
the data on the postglacial upheaval of Fenno- 
scandia. 


Our calculations were in the first place 
concerned with the question regarding the 
possibility of a relatively thin crust. We may 
now ask the question what possibilities for the 
explanation of the geological facts remain. 

First of all it might be assumed that the 
crust is thicker, i.e. the characteristic length 
for an elastic plate l is enlarged. As was shown 
in Part II, a value of l of 88 km, corresponding 
to a erust’s thickness of about 50 km was still 
not sufficient to explain the great distances 
between the positive zones and the Vening 
Meinesz zone. 

So we have to suppose this length to be 
some hundreds of kilometers. The gravity field 
in West Indonesia could then be accounted 
for by assuming a relatively light body in the 
crust, i.e. in the upper part of. it. For, 'by 
enlarging |, the thickness of the crust increases 
as well, and the cause of the negative anoma- 
lies cannot be situated very deep as a conse- 
quence of the steepness of the gradient. The 
positive anomalies would now be sufficiently 
far removed from each other to coincide with 
the positive zones on either side of the 
_ Vening Meinesz zone. 

A major objection against this construction 
however, is formed by the current estimates 
of the thickness of the elastic crust. A further 
explanation is now necessary, namely what 
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exactly must be understood by this term. For 
the elastic crust does not coincide with the 
seismic crust, which is clearly demonstrated 
by Ewing’s recent data", personaly commu- 
nicated to me by Professor Vening Mei- 
nesz, from which it appears that the Moho 
(rovicic) discontinuity, underlying the conti- 
nents at about 30 km depth (Veldkamp, 1949), 
lies at about 5 km depth beneath the ocean 
floor. With such a trifling thickness of the 
elastic crust no volcanic island would continue 
to exist longer than some tens of thousands 
of years. 

Nor need the isostatic compensation level 
be in direct relation with the elastic crust. 

The gravimetry, however, has still more 
sources of information at its disposal. From 
the application of regional isostatic reduction 
follows a value of l (see formula 5), the 
characteristic length of a floating elastic plate. 
From the formula it follows that the thickness 
of the elastic crust is determinative for this 
length. Vening Meinesz (194la) on the ground 
of comparative calculations about the degree 
of regionality at which the Hawaiian Archipe- 
lago is compensated, arrived at a value for | 
from 60 to 80 km and for the crust’s thickness 
T from 30 to 40 km. 


independently of Vening Meinesz, Ross 
Gunn (1943a and b, 1944, 1949) concluded 
in a series of articles, among which there is 
also an elaboration of the isostatic compen- 
sation of Hawaii (1943b), that b, a quantity 
that is to be put equal to (1.3 x 1)-!, attains 
a value of 84 x 10-8cm-1, which is about 
equal to 1 = 100 km. From this Gunn makes 
the conclusion that there is an elastic crust 
of 50 #15 km. 


Further evidence can be found in the 
structure of the African graben. Vening 
Meinesz (1950 b) in dealing with the pro- 
blem of this graben with regard to the question 
whether these structures have to be consi- 
dered as compression or tension phenomena 
— presuming that the elastice crust has a thick- 
ness of 40 km — stated that the calculated 
deformation is in good agreement with the 
actual width of the graben. !! This implies that 
the presumed thickness of the elastic crust 
(40 km) indeed is very likely. 


These data clearly show the untenability of 


ı0 This data will be published in the February 
issue of the Bull. Geol. Soc. America, Vol. 65, 1954. 

11 In this Vening Meinesz started from the data 
and calculations of Sir E. C. Bullard: “Gravity 
measurements in East Africa” (Phil. Trans. R.S., 
757, Vol. 235, pp. 445—531, 1936). 
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the supposition of l having so great a value. 
Wide perspectives are opened if we assume 
considerable horizontal forces working in the 
cerust which is supposed to have a thickness 
as stated above. These forces could easily 
account for the negative as well as the positive 
anomalies with respect to the regional positive 
field of about +20 mgals in Indonesia (see 
e.g. Vening Meinesz, 1940, 1948, 1950a). 
Moreover these forces would have a streng- 
thening effect on the earth’s material. Starting 
from the maximum shear theory it can be 
computed by means of formula (7), that at the 
surface a horizontal stress of 500 kg/cm? may 
occur, at a depth of 25 km a stress of more 
than 2500 kg/cm?. On the plane AB presented 
in fig. 11 these forces exercise a strengthening 
influence that is comparable with an increase 
of the hydrostatic pressure. 


It has been said before that the maximum 
shear theory is no sound basis for strength 
calculations. So the numerical values in the 
above cannot be nothing more than a guess. 


It may be pointed out here that Bijlaard 
(1933) as well as Vening Meinesz (1950a) 
later on presumed that a horizontal pressure 
of some thousands of kg/cm2 is needed for a 
buckling of the crust. The latter estimates 
these forces at about 2000 kg/cm2. 


Considering now the results of our compu- 
tations on deformations of the crust caused 
by essentially vertical loads, it may be 
concluded that the solution of the problems 
on mountain building must probably be sought 
for in assuming considerable primary hori- 
zontal pressures in the earth’s crust. Much 
work has already been done in this field; the 
above mentioned publications of Vening Mei- 
nesz and Bijlaard are some examples of this. 


It is beyond the scope of this study to give 
a detailed synopsis of the various hypotheses 
which are built on the assumption of hori- 
zontal pressure in the crust. It was our object 
to derive some boundary conditions for hypo- 
theses on mountain building in general. 

A first limitation of the possibilities has 
been given, We are fully aware that the last 
words on this question have not yet been 
spoken. Some, but certainly not all geophysical 
aspects of Indonesia were taken into considera- 
tion. Moreover, many other regions are open 
to similar investigations. 


Finally, some remarks must be made on the 
character of our computations. We usually 
started from schematic simplifications, for 
instance the earth’s crust was regarded as an 


elastic, homogeneous plate floating in a fluid 
substratum; a sharply defined transition from 
the crust to the substratum was assumed 
instead of a continuous one. Such simplifica- 
tions naturally affect our results. We are of 
opinion, however, that none of these simplifi- 
cations had an essential influence on our 
results and that abolishing them will not 
change the core of our argument. 

At any rate there seems to be no justification 
in building on possibilities that are at variance 
with the results of this study, without having 
first proved their physical acceptibility. 
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DEARAUR IE AINISEI PA EHE3IONTSN0:9, U N I9NSTA 


ON ONE OF THE OLDEST KNOWN REMAINS OF THE COMMON ELK, 
ALCES ALCES L., FOUND RECENTLY IN THE NETHERLANDS 


D. P. ERDBRINK 


SITUATION AND DESCRIPTION 
OF THE FIND 


In September 1953 the leaseholder of the 
loam pits near Hattem (province of Gelder- 
land, Netherlands), Mr. Keijl, sent me a notice 
that some bone fragments had been discovered 
during digging operations at a new spot in 
the old loam pit, designated as loam pit A 
in an earlier publication by De Waard and the 
present writer (1949). Upon investigation these 
"bone fragments proved to be still partly in situ 
" in a layer of ferruginous concretions at an 
approximate depth of 8 metres below the 
present surface near the pit. Unfortunately 
the find had been badly shattered because a 
lorry had driven over it. The bone splinters 
were lifted out carefully, as they appeared to 
be very soft and friable. 47 fragments were 
collected, cleaned and subjected to an immer- 
sion in a shellac solution, after which they 
could be fitted together in their” original 
position. 

The total absence of a brow tine made it 
apparent that this fossil fragment of an antler 
could only have belonged to an elk. The 
reconstruction of the object showed it to be 
the basal part (beam and beginning of the 
 palmated area) of the right antler. The con- 
tinuation of a small part of the bone below 
the burr proves that the antler was not a 
shed one but that it must have been part of 
a complete skull at the time of the animal’s 
death. The burr is not distinctly visible around 
" t+he whole circumference of the base of the 
beam, this being only the case around its 
lower half (or under side), but its original 
position can be inferred with sufficient cer- 
tainty so that one may decide that ır ıszthe 
right, and not the left, antler. The burr has 
evidently been subjected to corrosion, a thing 
which seems to happen easily in Alces, at least 
according to Guenther (1951, p. T20)rs, his 
corrosion can either have been a chemical one 
_ (due to being imbedded in a mass of ferru- 
ginous concretions), or one of a mechanical 
“nature, for instance caused by prolonged 
transport in water. During the cleaning of 
the lower part of the beam some curious 


depressions and striae, filled up with loam, 
became visible. The only explanation for these 
phenomena, in my opinion, is the assumption 
that they have been caused by scavengers, 
thus strengthening the supposition which I 
expressed earlier (De Waard and Erdbrink, 
1949, p. 686). 


From above, and to a lesser extent from 
below, the beam shows fairly strong parallel 
furrows, giving the whole a cannelated or 
fluted appearance. This same phenomenon is 
observable, on occasion, in the antlers of recent 
Alces; probably it is related with a fairly 
strong development of the blood arteries which 
have fed the antler at the time of its growth. 


The form of the palmated area of the antler 
cannot be inferred from the little which has 
been left. However, I think that it is fairly 
certain that it was not an antler of the so- 
called “cervine” type, such as has been defined 
by Lönnberg (1902) in a very important paper 
on the variability of elks’ antlers. In the 
Hattem fossil there have probably been at 
least four tines, while the type of antler must 
have been the palmated one or the one inter- 
mediate between palmated and cervine. 


The textfigures given by Lönnberg and 
others very clearly show that the form of the 
beam is as much subject to variability as the 
form of the antler, it being either straight 
or slightly to strongly curved downward and 
then again upward into the “palm” or into the 
point where the beam starts to branch. I have 
not been able to find many allusions to this 
feature among the authors whom I consulted. 
It seems certain, though, that Guenther’s 
assumption (1951, p. 122) that the beam is 
curved in most cases in Alces, does not appear 
to be in accordance with the facts which I 
collected from textfigures in literature and 
from comparison with a number of recent 
elk-antlers. 


Table 1, in which measurements of the fossil 
fragment from Hattem are compared with 
those of several recent and fossil antlers 
(mostly from various European localities), 
shows that here also there is nothing which 
could serve to distinguish the Hattem find 
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from the recent Alces alces !. I compared three 
measurements, to wit: 

(1) the circumference of the burr; (2) the 
minimum circumference of the beam; (3) the 
distance from the burr to the point where a 
broadening of the beam (the starting point 
for the palmation or the branching) becomes 
apparent, measured along a horizontal line 
on the upper side of the antler, and at right 
angles to the burr. 

Bächler (1911, pp. 81—82) has given a good 
description of the formation of the antlers in 
Alces, from which it follows that the rule of 
“one pair of tines added each year” is by no 
means true. Usually the antlers with two or 
more tines are only formed after the animal 
is at least four years old (they become adult 
after approximately five years). Old animals 
tend to develop antlers which are less imposing 
than those of specimens at the height of their 
existence, which is reached at an age of 16 
to 20 years. 


COMPARISON OF THE FIND WITH 
OTHER FOSSIL ELKS 

As has been remarked by Frentzen (1928, 
p. 385), there is no reason to suppose that 
the early members of the species alces were 
in any way more robust, or that they had 
larger antlers, than their recent counterparts. 
The fossil from Hattem, in no respect a large 
specimen, is in accordance with this remark. 
One could even suppose with Dietrich (cited 
by Frentzen) that the recent elks possess the 


1 As appears to be almost usual in Zoology, 
many synonyms of the name of the recent elk, 
Alces alces (L.), are in existence. One may en- 
counter, after Linne’s Cervus alces and Alce alces 
(a Latin misnomer; the singular of the word, as 
used by Pliny, is aleis, plur. alces): Alces palmatus 
Blasius aut Owen; Alces jubates Ogilby; Alces 
machlis Ogilby; Alces resupinatus Rouillier; Alces 
Savinus Rouillier; Alces antiquorum Rüpp; Alces 
bedfordie Lydekker; Alces diluvii Pohlig; etc. 
As a survey of the recent elk is not the object 
of this paper, I feel that I am not in a position 
to decide whether there remains any justification 
for keeping apart, as a separate species, the 
American elk (or Moose), under the name Alces 
americana, as is usually done. Those few spe- 
cimens and data wich came to my attention 
seemed to me to indicate the absence of distinguish- 
ing characters of any constancey between A. 
alces and A. americana (in view of the great 
variability of the recent elk in size, form of 
antlers, and in other respects). Moreover, the elk is 
an animal which often travels very long distances 
(see Bächler, 1911, in his description of the 
habits of the animal), a habit which is unfavour- 
able to the developing of separate races, sub- 
species or species. A more thorough study would 
probably reveal the untenability of A. americana 
as a Species. 


larger antlers when compared? with those 
of their species known from the Pleistocene. 
Frentzen mentions a set of antlers (coll. Graf 
von Ehrbach) with 41 tines as one of the 
largest known recent European elks; a set 
with 43 tines from the Kenai peninsula, Alaska, 
in the Duke of Westminster’s collection being 
a particularly large specimen of the North 
American elk (or “Moose”). Bächler (1911, 
p. 85) mentions the enormous antlers of an 
American elk in the Chicago Natural History 
Museum (then: Field Columbian Museum). 
From the deposits of the river Theiss (or Tisa) 
near Fegyernek in Hungary a fossil pair of 
elk’s antlers with more than 40 tines is known 
(National Museum, Budapest). 

A look at the table which accompanies this 
paper suffices to show that the fossil from 
Hattem is probably not an Alces latifrons 
R. Johnson, as the chief distinguishing charac- 
teristics of this early Pleistocene species seem 
to be the much bigger distance between the 
burr and the point of broadening of the beam, 
and the larger circumference of the beam®. 


2 The few data compiled by me in the table of 
measurements tend to support this view. One 
has only to compare the computed maxima and 
minima of the measured recent and fossil 
specimens of Alces alces. 

3 Schmidtgen (1938, pp. 334-335) described a 
curious abnormality, encountered by him in a 
shed left antler of A. latifrons from Mosbach, 
Germany. This antler possessed a kind of brow 
tine, starting at the burr. Apparently it (i.e, the 
tine) resembled, in itself, the form of a two- 
tined antler of a young (subadult) elk. Schmidtgen 
thought that this phenomenon has perhaps been 
caused by an injury to the burr of the animal, 
before its antler had started to grow, He had not 
been able to find anything of the same nature 
among recent elks. 

In this respect I have perhaps been luckier, 
because I encountered a nearly identical for- 
mation in a recent, stuffed specimen in the 
collection of the Rijksmuseum van Natuurlijke 
Historie at Leiden. This animal, according to the 
attached label, came from Canada and had lived 
in the Zoological Garden at Rotterdam for several 
years. It possesses 11 tines (5 in the right, 6 in 
the left antler), of which one tine in the left 
antler is curved downwards in a pathological 
manner. At a point 4 centimetres from the burr 
of the right antler, a tine juts out to the front, 
giving the impression of being a brow tine. The 
burr of the left antler has a short tine jutting 
out of it to the front, preeisely in the way 
described by Schmidtgen for the A. latifrons from 
Mosbach, and resembling the brow tine of a deer. 
The type to which the antlers of this recent 
specimen at Leiden belong, is the intermediate 
one between “cervine” and palmated. The state- 
ment on the attached label, that the animal was 
only three years old at its death, must be an error 
because of the way of development of the antlers, 
and the number of tines. 


Schmidt (1934) declared-“that these two 
measurements may be 1,9 times as much in 
"A. latifrons as in the recent elk (most skeletal 
| measurements of latifrons being about 1,5 
| times as large asin A. alces). It must, however, 
remain questionable whether this statement 
(based on finds made in Germany?) by 
Schmidt remains valid everywhere for A. 
| latifrons. The commonly accepted specifical 
distinctions of A. latifrons (compared with 


the recent elk) have been enumerated by . 


| Reynolds (1934, p. 2) as follows: (1) wider 
! frontals, (2) longer beam of the antlers, (3) 
4 simpler palmation of the antlers, (4) larger 
size in all respects. 


But at the same time he criticized these 
distinctions, stating that (in the Forest Bed 
specimens) “...in most of the specimens the 
beam is no thicker than in some recent 

“antlers.” The palmation of the specimens 
| known to Reynolds was not very great, and he 
concluded that ”...the most valid characters 
of Alces latifrons would seem to be the simple 
palmation and the length of beam.. Dabutisne 
remarked that there is no palmation in the 
“Eastern Siberian variety of EIk”. Here again 
we may revert to the opinion expressed by 
Lönnberg (1902, p. 353), who showed that the 
supposed specifical characters of some elks 
from Eastern Siberia, described by Lydekker 
(1902) as Alces bedfordie, were in reality 
nothing else but examples of the cervine type 
of antler, without palmation, which is also 
common in Scandinavia and the Baltic coun- 
tries. As a result of these considerations Ithink 
‚that the only apparent difference which serves 
to keep apart the antlers of the fossil A. 
latifrons, and the fossil and recent A. alces, 
seems to be the greater length of beam of the 
first species. 


It may be remembered here that some other 
differences between the two species are known, 
all of which tend to show that A. latifrons 
was not so well adapted to life in wooded and 
swampy country a8 the present A. alces. 
Schmidt (1934) observed that the bones of the 
feet (especially the phalanges, metacarpalia 
and metatarsalia) of latifrons indicate that 
this animal did not possess the broad hooves 
which serve the recent elk so well when 
traversing swamps and streams. 


I have not been able to find any data on 
the presence or absence, in latifrons, of the 
snout-like apparatus of nose and upper lip, so 
characteristical to the recent A. alces, and 
indicated by a shortening of the nasalia in 
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the skull. Thus far this shortening has not been 
described in any skull of A. latifrons. 

Lönnberg (1902, p. 360) concluded that “the 
forked stage of the immature antler with 
rounded tines may to a certain extent be 
regarded as a repetition of the phylogenetic 
development, so that in this way the “cervine” 
elk-antler, whether it be called a development 
of the young stage or a reversion, displays 
primitive characteristics in its rounded tines.” 
The simpler palmation, usual in A, latifrons 
(although the specimen figured by Schmidt, 
1934, has magnificently palmated antlers), 
seems to be consistent with Lönnberg’s remark. 

As was pointed out by Schmidt, Soergel 
(1914), as a result of a revision of the faunae 
from a number of middle Pleistocene German 
localities, concluded that there appears to be 
a well-marked, but slight, difference between 
the specimens of A. latifrons from Mosbach 
and those from Mauer. The latter ones are 
smaller in most respects, some falling within 
the range of measurements of the recent elk. 
Soergel concluded that the typical A. latifrons, 
such as those from Mosbach, have been 
inhabitants of the steppes*, whereas the 
animals known from Mauer have been living 
in a wooded biotope. He thought that these 
smaller elks, inhabiting woods as does the 
recent elk, must be placed in the direct line 
of ancestors of A. alces; A. latifrons (type), 
the steppe animal, having died out without 
leaving any direct descendants. 


STRATIGRAPHICAL OCCURRENCE 
OF THE GENUS ALCES 


The stratigraphical location of the last 
members of the species latifrons is commonly 
placed at the beginning of the Penultimate (or 
Mindel-Riss) interglacial, as is apparent from 
the descriptions and table given by Hescheler 
and Kuhn (1949, p. 175). The deposits at the 
German locality of Süssenborn (dating from 
the immediately preceding Antepenultimate 
or Mindel glaciation, see Erdbrink, 1953, fig. 
37, p. 432) still contain many remains of 
A.latifrons, according to Soergel (1914, p. 148), 
Schmidt (1934), and others. 

The first occurrence of specimens of the 
recent elk, A. alces, according to Hescheler 
and Kuhn (1949, p. 173), Guenther (1951, p. 
120), Bächler (1911, p. 100) and others, is 


4 Roman and Dareste de la Chavanne (1931, 
pp. 1256—1257) reached this same conclusion with 
regard to an Alces latifrons from Seneze (dep. 
Haute-Loire, France), of which they superficially 
described a nearly complete and very large 
skeleton. 


ran een 
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nr 5 TABLE OF MEASUREMENTS 
(in centimetres) 


Distance nina 
burr-be- mMMUM| Orcum- 
Recent specimens: ginning of en f ference | Remarks 
palm or ee of burr 
fork 
Rijksmuseum v. Nat. Hist,, r: 89 14,2 22 Palmated antlers; 14 tines. Vertical 
Leiden cat. IX, a. I 95 13 21,3 section of beam circular. Lithuania. 
Rijksmuseum v. Nat. Hist., r 10,7 12,8 17,5 Cervine antlers; 6 tines. Vert. sect. 
Leiden cat. IX,, d. ] 11,3 12,9 17,9 of beam circular. Lithuania. 
;: : x i Cervine antlers; 4 tines; degenera- 
Rijksm ; h h & 4 Ä 
A end 1 | ns ee os | ted. Vert. sect. of beam oval. Lithu- 
ÄcE i : 3 i ania. 
Rijksmuseum v. Nat. Hist., T“ 95 14,2 21 Palmated antlers; 11 tines. Coll. 
Leiden cat. IX, £. 1. 9,5 14,7 22.5 Joh. Paludanus. 
Rijksmuseum v. Nat. Hist., D: 16 18 30,5 Palmated antlers; 33 tines. Vert. 
Leiden cat. IX, g. % 18,5 lin 29,2 sect. of beam circular. Lithuania. 
Rijksmuseum v. Nat. Hist., T, 12 18 27 Palmated antlers; 21 tines. Vert. 
Leiden cat. IX, h. 1. 12,5 18,8 26,5 sect. of beam circular. Lithuania. 
Rijksmuseum v. Nat. Hist., ig 11 15,3 23,2 Intermediate type. 14: tines. Vert. 
Leiden cat. XI, a. j. 18 16 95 sect. of beam circular. Coll. Thun- 
berg, Lapland. 
Rijksmuseum v.Nat.Hist,| r. 4 17,2 Bo laer medialen I u Yen 
een CMIXT IE: 1 9 164 22.5 sect. of beam oval. Pathological case 
g ; : “browtines”). Canada. 
« : Cervine antler, shed (Zoo, Amster- 
an e: as = ze | u. z_ dam; don. Dr. A.E. v. Giffen). Vert. 
2 . ; sect. of beam circular. 
eg v Arte 1. u 18 30 Two shed left antlers, Lithuania. 
ELF EN Are a. has 7 tines and circular vert. 
Id., reg. no. 2548, b. E 11,5 19,5 30,3 sect. of beam, b. 8 tines and oval 
sect. of beam. 
Coll. Miss I.M. Saetrang, r 11,6 16,6 21,7 
Baarn 1 18,7 16,6 22,6 Palmated antlers. Vert. sect. of 
: beam oval. Both features identical 
Coll. J. Smit, Baarn (see T. 10 15,8 221 with 2nd spec. also from Sweden. 
phot.) il 319 15,6 21,4 
Zool. Inst., Kiel, 16382 T} 16 16 26 
(Guenther, 1951) 1 15 16 25 
Geol. Inst, Kiel (Guen- r 13 13,5 20 
ther, 1951) il 13 14,5 23 
Karlsruher Landessamm- T. 4 — — Odd 14 tine-specimen from Norway. 
lung (Frentzen, 1928) 1 5 — — 
Karlsruher Lande San. 7; 9,3 _ _ 16 tine-specimen from Norway. 
lung (Frentzen, 1928) Naz 8,8 = = 
Karlsruher Landessamm- r 123 — = Odd 18 tine-specimen from British 
lung (Frentzen, 1928) 1 14 — — Columbia. 
Karlsruher Landessamm- > 13 —_ — 28 tine-specimen from British 
lung (Frentzen, 1928) 1: 13 — u Columbia. 
Fig. 1 — (a) Fossil fragment of an antler of Alces alces L., seen from above. Found at a depth of 


8 metres below the surface in loampit A (the „Oude Leemkuil”) at Hattem, prov. Gelderland. 
(b) The same fragment, seen from below; remains of the burr are visible. 
_ (e) Part of an antler of a recent elk from S. Sweden (coll. J. Smit, Baarn), seen from above; 
compare with fig. la. 
(d) Part of the same recent, right antler, seen from below; compare with fig. 1b. 
All figures approximately % nat. size. 
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TABLE OF MEASUREMENTS 
(in centimetres) 


Distance le | 

burr-be- eng Circum- 
Fossil and subfossil SUR eircum- ference Remarks 

i ginning of ference of 
specimens: palm or ee of burr 
fork 
Zool. Inst., Kiel, sub 62 1. 115 15 20 Site unknown; postglacial. 
(Guenther, 1951) 
Zool. Inst., Kiel, sub 74 I ital 15,2 24 Torfmoor Stormarn; subrecent (?). 
(Guenther, 1951) 
Zool. Inst. Kiel (Guen- | 2 14 20 27 Moorkirchen, Westerholz; 
ther, 1951) ik 12 20 27 postglacial. 
Zool. Inst., Kiel, sub 75 1, u 20 20 | Warleberger Moor, postglacial. 
(Guenther, 1951) ll, 10 19 j 27 
Geol. Inst, Kiel (Guen- T% 18 18,5 = 769 Preetz; Weichsel interstadial. 
ther, 1951) il 16 18,5 =22 
(Frentzen, 1928) Ta 118} 14 16,4 Lower Rhine terrace, Knielingen; 
: young Pleistocene. 

& Rügetswil-Junkertswil, Switzer- 
at) 2 15,8 19,6 a land; young Pleistocene or post- 
il 14,8 12,6 R at 

Geol. Inst., Utrecht, Ik 12,5 at 20,9 E, Scheldt; young Pleistocene. 

G 21—1943 

(De Man, 1878) 1 5 — - W. Scheldt near Breskens; young 
Pleistocene. 

Coll. J. Butter, Deventer m 15,5 12 E 22,5 New harbour along river Yssel, 
Deventer; young Pleistocene. 

Maxima & Minima 18 —10 120 —12 127 —164| 

(fossil specimens) 

Maxima & Minima 18,5— 4 119,5—10,5 i30.5—12,5 

(recent specimens) 

The described specimen S Et ESE = 12,8 

from Hattem 

Alces latifrons ik 24 22,5 36 From Mosbach n. Mainz (abnormal 

(Schmidtgen, 1938) antler), 

Alces latifrons 2 46 — m From Mosbach n. Mainz (normal 

(Schmidtgen, 1938) antler), 


usually given as the beginning of the Last 
(or Riss-Würm) interglacial. I have not been 
able to find any allusions to older finds of 
A. alces in the literature consulted by me. If 
one accepts Soergel’s theory about the pri- 
mitive elk from Mauer (upper part of the 
Antepenultimate or Günz-Mindel interglacial) 
being the direct ancester of A. alces, there 
exists a time lapse from the Antepenultimate 
glacial till the Last interglacial during which 
nothing is known of A. alces or its direct 
ancestors. Even when the — extremely impro- 
bable — theory is accepted that the recent elk 
has had A. latifrons as its direct ancestor, no 
remains of the genus Alces have been descri- 
bed which can be dated back to the Penultimate 
(Mindel-Riss) interglacial, or to the Penulti- 


mate (Riss) glaciation; however, there is one 
exception, namely the mention made by 
Oakley (1938, p. 28) of remains of an Alces 
sp. (tentative identification) occurring in the 
Middle Gravels of the Barnfield Pit at Swans- 
combe in the Thames valley, a deposit which 
has been dated to the Penultimate (Mindel- 
Riss) interglacial. Elsewhere in Europe the 
absence of Alces alces prior to the Last inter- 
glacial is apparent, as is testified by many 
publications about finds of fossil remains of the 
recent elk (Bächler, 1911; Degerbgl, 1933, pP. 
359, 365, 372, 374, 385; Frentzen, 1928, pp. 381— 
385; Guenther, 1951; Reynolds, 1934; Rüeger, 
1944, pp. 282—283). 

In early historical times A. alces has had a 
much wider distribution than is commonly 


known. According to data'mentioned by De 
Man (1878), and by Bächler, in a very thorough 
survey of the fossil elks found in Switzerland 
(1911; pp. 96 et seqg. contain a lot of useful 
information about the habits of the recent 
elk), A. alces was still an inhabitant of the 
Netherlands and northern Belgium in the 
tenth century _A.D.>. 


OTHER FOSSIL ELKS FROM THE 
NETHERLANDS 


- Fossil remains of Alces alces have been 
“found in the Netherlands, but it seems that 
all of these, with reasonable certainty, can be 
dated back to more recent times. De Man 
(1878) gave a detailed description of a left 
antler of an elk, dredged at a depth of 20 to 
30 metres from a narrow channel between two 
mudbanks along the coast near Breskens 
.(opposite Flushing) in the Western Scheldt. In 
view of the other finds which he described 
from the same area (among which are abun- 
dant remains of the Mammoth) and of what 
is now known about the fossils from the 
Scheldt estuary, we may safely assume that 
this elk belonged to a fauna living near those 
regions during the Last (Riss-Würm) inter- 
glacial or the Last (Würm) glaciation. The 
very fine engraving of the find, published by 
De Man, shows an antler of the intermediate 
palmated-cervine type, with 5 tines. This 
refutes Bächler’s (1911, p. 84) statement that 
fossil antlers of the cervine (or the inter- 
mediate; “Stanglergeweihe”) type have never 
been found ®, 

Rutten (1909, pp. 68—70) mentioned and 
shortly described four other fossil or subfossil 
remains of elks from the Netherlands. One, 
a nearly complete skull in the collection of the 
Rijksmuseum van Geologie at Leiden (St. 
12288), and found near Lisse, had already been 


5 As is testified by a charter from the year 943, 
given by the Emperor Otto I (the Great), 
mentioning the exclusive rights of bishop Balderic 
of Utrecht to hunt elks in the forests of Drenthe 
(a province in the northern part of the Nether- 
lands). A Flemish chronicle of the tenth century 
mentions the elk (often called the elo or schelo, 
which is an ancient Saxon word; compare Sieg- 
fried’s hunting prowesses in the Nibelungenlied, 
St. Gallen manuscript, stanza 3753, according to 
Bächler: elch, Schelch) as an inhabitant of 
Flanders. 


6 Among the many excellent photographs ac- 
companying Bächler’s paper, there is one (Pl. VII, 
fig, 2) of a fossil right antler of Alces, found 
at Gossau, Switzerland, and which is of the 
cervine type. 

One of the two specimens found at Deventer, 
and in the collection of Mr. Butter (number 
D28 NH), is of the intermediate type. 
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mentioned by De Man (1878). These remains 
all came from early Holocene or late Pleisto- 
cene deposits, certainly not older than the 
Last interglacial. Kunst (1937, p. 29) added 
one other find to those mentioned by Rutten; 
it came from Hengelo and was found in the 
course of digging operations for a lock in a 
canal. It is of Last (Würm) glacial age. 


In addition to these finds, the Mineralogical- 
Geological Institute of the Utrecht University 


“ possesses a fossil left (shed) antler (No, G 21, 


1943; don. Dr. P. Vermey), dredged from the 
Eastern Scheldt near Zierikzee, and also 
datable to a probable Last (Würm) glacial 
age. A set of fossil antlers of Alces alces, found 
during dredging operations in the river Yssel 
near Zutphen, is in the townhall at Zutphen 
(oral communication by Mr. J. Butter); its 
age, and that of two other specimens in Mr. 
Butter’s collection at Deventer (a find made 
in the course of the construction of the new 
harbour and a lock giving access to it, along 
the river Yssel at Deventer) is certainly not 
older than the Last interglacial, and very 
probably younger. 


PARTICULARS OF THE FOSSILIFEROUS 
DEPOSIT AT HATTEM 


The fossil antler fragment from Hattem, 
which should be identified as a specimen of 
Alces alces (L.), and which has been added 
to the collection of the Rijksmuseum van 
Geologie at Leiden under the number St. 28379, 
is of some importance because it helps to fill 
the stratigraphical and phylogenetical gap in 
the history of the genus Alces. The loam 
deposits at Hattem have been deposited during 
the initial part of the Penultimate (or Riss) 
glaciation, according to De Waard and Erd- 
brink (1949). Contrary to my statement in 
that paper, that the six fossil fragments of 
antlers (among which there is one piece, c, 
which might also be ascribed to an Alces sp., 
as was remarked by Hooyer) were found 
directly below the varved series, thus im- 
plicating the absence of varved loam lower 
down, it has now become clear that these 
varved series have a larger downward exten- 
sion in loampit A. The “lowest” layer of loam 
with the peculiar blue-grey colour, followed 
(in a downward direction) by the layer of 
ferruginous concretions which contained the 
fossils, and under which is a layer of whitish 
sand of much greater thickness than the other 
layers in the pit, were met with again in the 
newly dug trench in the bottom of loampit A. 
But because this trench has been dug down, 
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at an angle, to a considerable depth, it now 
appears that below the layer of white sand 
(which may reach a thickness of 50 centi- 
metres) the varved series start anew. I counted 
at least 70 separate thin layers of loam the 
last time when I visited the pit (December 4#, 
1953), while digging operations in the inclined 
trench were still going on owing to the 
presence of yet more loam lower down. The 
comparatively thick layer of white sand,. and 
the layer of ferruginous concretions must 
therefore have been deposited, not at the start 
of the deposition of the varved loam (indi- 
cating the nearby presence of an icecap of 
the Riss glaciation), but during an interruption 
of this deposition, which was perhaps due to 
some local or temporary regression of the 
icecap further to the North. This would also 
explain the presence of these animal remains 
in an otherwise completely sterile deposit of 
glacial origin. It is very peculiar that up till 
now only fragments of antlers have been 
recovered from the ferruginous layer, and no 
other bones or teeth. The presence of some 
kind of selective influence at the time of their 
deposition (scavengers, man?) cannot be 
completely ruled out. A curious fact in relation 
to the occurrence together of antlers of the 
red deer, Cervus elaphus, and Alces alces at 
Hattem, is formed by Bächler’s statement 
(1911, p. 86) that these two species seem to 
detest each other, avoiding each other’s 
presence as much as possible. 


At the spot where the badly shattered antler 
fragments of Alces were found, I measured the 
strike as being N 30°W, the dip (to the West) 
being 27 degrees at the same place; but at 
about 30 feet straight to the South, the same 
(ferruginous) layer had a strike of N 84°E 
and a dip (to the East) of 22 degrees! This 
shows how the whole deposit of loam and 
sand has been kneaded and broken into 
separate miniature fault blocks. Independent 
movement of the sand and the loam layers 
may perhaps be held to be one of the abrasive 
influences responsible for particular spots on 
some of the fossil antlers, where there are a 
number of peculiar flattened regions (De 
Waard and Erdbrink, 1949, phot. 9, a, e; and a 
hitherto undescribed new fragment of an 
antler, found by me recently at Hattem, and 
perhaps identifiable with C. elaphus). 


Concluding, it may be maintained that, in 
the absence of more definite data about the 
Alces sp, from Swanscombe, the antler frag- 
ment from Hattem must be seen as one of the 


oldest known remains thus far of the common 
elk, Alces alces, from anywhere in the world. 
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»DEFTEERBOGCELITZ 


NIEUW ZEEFAPPARAAT VOOR HET NAT AFZEVEN OP KLEINE 
ZEEFMATEN 


F. J. FONTEIN ! EN H. DREISSEN 1 


INLEIDING 


* Op het Centraal Proefstation van de Staats- 
mijnen in Limburg werd kort geleden een 
apparaat ontwikkeld voor het nat afzeven op 
kleine zeefmaten. Dit apparaat, zeefbocht ge- 
naamd, wordt in de wasserijen van de Staats- 
mijnen reeds op diverse plaatsen toegepast 
voor het afzeven en ontslikken op 0,5 mm. 


| VOEDING 


TOEVOERPYP 


VALTRECHTER 


VERDEELBAK 


500 


SPLEETZEEF 
DAMBREEDTE 
CA 2,0 mm 
SPLEETWUDTE 
CA1,0mm 


DOORVAL 


Fig. 1 — Zeefbocht voor afzeving op 0,5 mm 
(maten in mm). 


BESCHRIJVING 


Een zeefbocht geschikt voor afzeving op 
0,5 mm vindt men op fig. 1 in doorsnede ge- 
schetst. Het af te zeven mengsel van deeltjes 
en water wordt in de verdeelbak toegevoerd 


1 Wasserijtechnische Sector, Centraal Proef- 
station van de Staatsmijnen in Limburg (Tree- 
beek). 


en stroomt via de overloop van deze ver- 
deelbak en de valtrechter tangentieel op 
een gebogen zeefoppervlak, dat bestaat uit 
roestvrij-stalen spleetzeef met een spleet- 
wijdte van ca. 1 mm en een dambreedte van 
ca. 2 mm. De ca. 2 mm brede profielstaven 
van de spleetzeef hebben in verband met de 
slijitage van de zeef bij voorkeur een recht- 
hoekige doorsnede en staan loodrecht op de 
richting waarin de vloeistof met de vaste deel- 
tjes over de zeef stroomt. 


VERKLARING VAN DE WERKING 


De zich snel over de zeef bewegende laag 
van water en deeltjes wordt bij het passeren 
van iedere spleet van 1 mm ca. 0,25 mm 
dunner, daar een laagje van deze dikte door 
de spleet naar buiten spuit. Door dit laagje 
kunnen deeltjes van 0,5 mm nog worden mee- 
gesleurd. Grovere delen gaan over de spleet 
heen, daar ze zich voor het grootste gedeelte 
in de vloeistof bevinden, welke over de spleet 
heen gaat. De spleetzeef is gebogen om er voor 
te zorgen, dat de vloeistoflaag met het zeef- 
oppervlak in contact blijft. Het is uit het voor- 
gaande duidelijk, dat door de spleten van 
l mm wijdte alleen deeltjes uittreden met een 
diameter kleiner dan 0,5 mm. Het apparaat 
zeeft dus af op ca. 0,5 mm en de spleten van 
l mm kunnen niet verstoppen, daar slechts 
delen van 0,5 mm de spleten passeren. Daar 
verstopping is uitgesloten, is de capaciteit van 
een zeefbocht onafhankelijk van de concen- 
tratie aan vaste stof en de korrelgrootte van 
de af te zeven delen. Alleen het aantal te 
verwerken m? per tijdseenheid is van belang. 


CAPACITEIT 


Daar de snelheid van de laagjes vloeistof, 
die door de spleten naar buiten treden, onge- 
veer gelijk is aan de snelheid, waarmede het 
mengsel van vaste deeltjes en vloeistof over 
de zeef spuit, is de capaciteit van een zeef- 
bocht recht evenredig met de toevoersnelheid. 
Bij een toevoersnelheid van ca. 3 meter per 


1, 9324 


seconde, optredende bij een lengte van de val- 
trechter van 50 cm, is de capaciteit van een 
zeefbocht volgens fig. 1 met een breedte van 
60 cm ruim 100 m3/h. 

De capaciteit per vierkante meter zeef- 
{ oppervlak van een snellopende schudzeef is 
bij ongeveer gelijke afzeving 10 tot 30 keer 
zo gering als die van een zeefbocht. 


ONGESLETEN ZEEF VOEDING 


DOORVAL 


DOORVAL 
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De spleetzeef moet bij toepassingen in kolen- 
wasserijen om de 1 a 2 weken worden ge- 
draaid. De constructie van de zeefbocht is dus- 
danig, dat de spleetzeef door het wegslaan van 
4 spieen los komt en gedraaid kan worden. Na 
ca. 4 maanden moet bij toepassing in de was- 
serij de spleetzeef worden vernieuwd, daar de 
afzeving te grof wordt. 


VOEDING GESLETEN ZEEF 


FUNERE AFZEVING 
MINDER GAPACITEIT 


OVERLOOP 


Fig. 2 — Invloed van de slijtage op de afzeving. 


SLIJTAGE 


Na enige tijd in bedrijf geweest te zijn, slijt 
de spleetzeef als op fig. 2 is aangegeven. De 
afzeving wordt dan vanzelfsprekend fijner. 
Door de spleetzeef te draaien, zodat de af- 
voerzijde toevoerzijde wordt, wordt de af- 
zeving weer grover. Door voldoende vaak de 
spleetzeef te draaien kan de fijnheid van af- 
zeving voldoende constant worden gehouden. 


GEGEVENS OVER DE RESULTATEN 
VERKREGEN MET DE ZEEFBOCHT 


Tabel 1 geeft als voorbeeld alle van belang 
ziinde gegevens betreffende een toepassing 
van de zeefbocht in een wasserij van de 
Staatsmijnen, waarbij de aftap van een spits- 
kast aan een zeefbocht werd toegevoerd. Zoals 
men zal opmerken is het toevoerproduct zeer 
moeilijk af te zeven, omdat veel toevoer- 


TABEL 1 — RESULTATEN VERKREGEN MET ZEEFBOCHT VOLGENS FIG. 1, BREEDTE 60 CM 


orrelverdeling 


Korrelgrootte Korrelverdeling RK 
u toevoer, % 


overloop ‚% 


Schets afmeting 
van profielstaven 


Korrelverdeling 
doorval, % 


2,0 


. 


< 2000 1,44 4,5 =; 
1000—2000 5,42 16,7 — Er 
710—1000 1urftt 5,8 0,1 \ 
500— 710 13,81 31,9 I) 
350— 500 16,43 19,7 14,9 
210-— 350 13,85 8,0 16,6 m 
105— 210 12,46 4,1 16,4 Maten in mm 
< 105 34,82 10,0 46,7 Afstand tussen 
profielstaven 
l mm 
Totaal 


Cap. in m3/h 
Conc. vaste 
stof in g/l 
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materiaal een korrelgrootte bezit, die slechts 
weinig afwijkt van de korrelgrootte waarop 
moet worden afgezeefd. 
CONCLUSIE 
De zeefbocht zal door haar grote eenvoud, 


VERSLAGEN VAN 


goede zeefeigenschappen en lage bedrijfskos- 
ten zeker vele toepassingen vinden in allerlei 
industrieen waar mengsels van vaste deeltjes 
en vloeistof worden verwerkt. Voor deze uit- 
vinding is octrooi aangevraagd. 


VERGADERINGEN 


VERSLAG VAN DE EXCURSIE NAAR DE TUNNELPUT VELSEN 


Op 22 Februari vond de door de Geologische 
Sectie georganiseerde excursie naar de bouwput 
voor de tunnelwerken bij Velsen plaats. 

In verband met de grote belangstelling, onge- 
veer 130 aanwezigen bij de inleiding en 100 bij 
de excursie zelf, werden de inleidingen tot de 
excursie, die om 10 uur des morgens in de biblio- 
theek van de Geologische Dienst geprojecteerd 
waren, nu gehouden in de grote gehoorzaal van 
Teyler’s Stichting, eveneens aan het Spaarne te 
Haarlem. Na de opening van deze dag door de 
voorzitter van de Geologische Sectie, Ir. A. van 
Weelden, werden een zevental korte inleidingen 
gehouden door Dr. J. D. de Jong (algemene 
inleiding), Ir. H. Wieringa van de Rijkswaterstaat 
te Velsen (ontwerp en uitvoering tunnelbouw 
Velsen) en voorts naar een stratigrafische in- 
deling van het profiel in de tunnelput door Prof. 
Dr. Mr. F. Florschütz (de pleistocene afzettingen), 
Ir. J. Bennema (bodemkundige aspecten van het 
laagterras), Drs. J. W. Chr. Doppert (het veen- 
op-grotere-diepte), Dr. L. M. J. U. van Straaten 
(de wad-sedimenten) en Ir. L. J. Pons (de jJongste 
afzettingen). Hierbij traden de sprekers mede 
namens de andere onderzoekers op, die tezamen 
deel uit maken van een werkgroep waarin de 
verschillende hier toepasseliijke facetten van 
geologisch, paleontologisch, bodemkundig en 
archaeologisch onderzoek zijn verenigd. Het list 
in de bedoeling de resultaten van al deze onder- 
zoekingen in &en grote Velsen-publicatie te ver- 
enigen. De lunch werd in het gebouw van de 
Geologische Dienst gebruikt. 

In Velsen werd het gezelschap in vijf groepen 
verdeeld, terwijl Prof, Florschütz en de heren 
Bennema, de Bakker (Stichting voor Bodem- 
kartering), Doppert, Van Straaten en Pons zich 
op verschillende plaatsen in de put hadden ge- 
posteerd en daar voor ieder der vijf groepen ter 
plaatse lieten zien wat des morgens reeds in een 
breder verband door hen naar voren was ge- 
bracht. 

Prof. Florschütz toonde het pleistocene deel 
van het profiel. Doordat het werk wegens vorst 
al vier weken stil lag en de vele wind in die 
periode in het zandprofiel de fijnere en grovere 
lagen had uitgeprepareerd waren hier zeer fraai 
de kris-kras gelaagde zanden met samenspoe- 
lingen van plantenresten, waarin Dryas voor- 
komt, en vorstspleten waar te nemen. Opmerke- 
lijk is het verschil met de horizontaal gelaagde 
zanden boven in het profiel; tussen beide zanden 
ligt een laag van kryoturbate kleien. 


Door de vorst was men met het werk nog niet 
zover gevorderd, dat de als Eemien aangeduide 
lagen ontsloten waren. 

De heer Bennema toonde in de kris-kras ge- 
laagde zanden onthoofde bodemprofielen. Deze 
vondsten waren nog van te recente datum om 
aan de hand hiervan tot een mogelijke uitspraak 
over een mogelijke Würmstratigrafie te geraken. 
In de top van het pleistocene pakket ziin bodem- 
profielen ontwikkeld, waarvan het karakter al 
naar de hoogteligging verschillend is. 

De heer De Bakker toonde in de oostwand van 
de put een niveau, dat sprekende gelijkenis ver- 
toont met het reeds in Noord-, Oost- en Zuid- 
Nederland herkende Allergd-niveau. Nadere aan- 
wijzingen ontbreken echter nog. 

De heer Doppert toonde het Veen-op-grotere- 
diepte. In dit veen komen klei-inschakelingen 
voor. Op dit veen, dat in de wanden van de put 
op de diepte waarop het verwacht kan worden 
overal wordt aangetroffen, rust een 1 a l% m 
dikke kleilaag, zeer rijk aan Hydrobia. Deze 
kleilaag vormt samen met het veen-op-grotere- 
diepte de scheiding tussen het brakke bovenwater 
en het zoete spanningswater. Het brakke boven- 
water wordt onmiddellijk gevoed door het Noord- 
zeekanaal; het zoete spanningswater maakt deel 
uit van de oostelijke uitbreiding van de zoet- 
waterzak onder de duinen. 

De door de heer Van Straaten te tonen geul- 
afzettingen in wadfacies bleken, in tegenstelling 
tot de andere ontsluitingen op deze dag, onhan- 
delbaar doordat zij nog een harde ijsmassa vorm- 
den. Het algemene karakter met de fijne klei- 
zand gelaagdheid was wel goed zichtbaar maar 
de fijne texturele verschillen om tot een facies- 
analyse van deze afzettingen te geraken konden 
nu niet gedemonstreerd worden. 

In de zuidwand van de put was een ontsluiting 
semaakt waar de heer Pons het contact van de 
geulafzettingen (fase Velsen II) met de oudere 
afzettingen (fase  Velsen I) kon tonen. 

Dankzij een archaeologische vondst kon de 
geul op enkele eeuwen vöör het begin van onze 
Jaartelling gedateerd worden. Tegen een vereen- 
zelviging van de kwelderklei met de Oude Zee- 
klei verzet zich een C 14-bepaling, verricht op 
schelpen (ouderdom ong. 3800 jaar) en de hoger 
dan normale lisging (top op 2.95—N.A.P.). De 
oudere afzettingen &n de geulopvulling worden 
afgedekt door zanden waarin veel Romeins is 
sevonden. Tenslotte wordt het profiel hier afge- 
dekt door lagen veenslik, zand en klei die voor- 
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lopig als post-16e eeuw gedateerd worden. De klei 
werd aan de oever van het IJ afgezet. 

Het voor West-Nederland gebruikelijke profiel: 
pleistoceen zand — Veen-op-grotere-diepte — 
Wadafzettingen — Oude zeeklei-duinzand of 
veen/jonge zeeklei is hier dus slechts in de eerste 
twee etages in zijn normale beeld ontsloten. De 
wadafzettingen en de oude zeeklei in de zin van 
de Geologische Kaart worden hier niet aange- 
troffen. Inplaats daarvan worden hier jongere 
afzettingen weliswaar in wadfacies aangetroffen, 
waarvan de oudste wellicht verbonden zullen 
kunnen worden met de opvulling van het Oer-IJ 
(Güray, 1951) en de jongere met een scherp 
contact een jongere inbraak-fase vertegenwoor- 
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Polarized light in metallography. Editors G. 
K. T. Conn and F. J. BrapsHmaw. ix + 130 
pp., 31 figs., 22 pls. Butterworths Scientific 
Publications: London, 1952, Prijs geb. 21 s. 


De toepassing van gepolariseerd licht bij het 
microscopisch onderzoek van metalen is betrekkelijk 
recent en nog niet algemeen, neemt evenwel aan 
belangrijkheid toe naarmate de ervaringen daar- 
mede groter worden. Het is daarom zeer te waar- 
deren dat met dit boek een zakelijk en goed over- 
zicht wordt gegeven over de tot dusver bereikte 
resultaten. 


De eerste twee hoofdstukken geven in beknopte 
vorm de theorie en de kennis daarvan is voor het 


digen. Wat in de wand van de bouwput boven de 
Hydrobia-kleilaag is ontsloten behoort dus be- 
halve voor de zuidwesthoek geheel tot de jongere 
geulafzettingen (enkele eeuwen vöör Chr.). Bij 
het profiel van 1941-’42 was reeds van deze 
afzettingen de westelijke helling in de Noord- 
wand door Doeglas als problematisch aangeduid. 
Om vijf uur werd de excursie te Haarlem ont- 
bonden. Het weer heeft buitengewoon medege- 
werkt om deze excursiedag te doen slagen. Was 
er daags te voren door een sneeuwdek weinig 
te zien van de vele interessante zaken, die het 
profiel van deze bouwput geeft, een milde regen- 
bui des nachts had de geologie letterlijk ontsloten. 
J. D. de Jong 
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verder doorwerken van het boek noodzakelijk. De 
apparatuur voor en het practisch onderzoek met 
gereflecteerd gepolariseerd licht worden in de hoofd- 
stukken 3 t/m 5 behandeld, waarbij duidelijk op de 
moeilijkheden en beperkingen gewezen wordt. Spe- 
ciaal hoofdstuk 5 over de insluitsels is voor het 
metaalonderzoek van eminent belang en zoals de 
schrijver terecht opmerkt ligt hier nog een groot 
terrein voor studie open. De juiste toepassing en 
interpretatie is moeilijk en vergt van de onderzoeker 
een grote mate van critische zin. 

Het boek kan de 
aanbevolen worden. 


serieuze onderzoeker warm 


A. F.E. Jomsen 
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De leden van het Genootschap worden verzocht 
alle wijzigingen in adres, titel en functie te willen 
opgeven aan het secretariaat: Mijnbouwstraat 20, 
Delft. 


Nieuwe leden: 


LINGEN, G. J. VAN DER — Maarssen, Binnen- 
weg 70 (bg) (U.G.V.) 


Bezwaren tegen de toelating moeten onder- 
tekend en met redenen omkleed binnen vier 
maanden worden ingezonden aan de Secretaris 
van het Genootschap, Mijnbouwstraat 20, Delft. 


Nieuwe Adressen: 
INDIEIVASE DIE GT 
straat 51 (g) 
BOERFINCHDENF Utrecht, Munststraat 3 (bg) 

(U.G.V.) 
BOISSEVAIN, Dr. H. — Wassenaar, Rijksstraat- 
weg 612 (g) (gk) 


— Amsterdam, Valerius- 


BOTTINGA, geol. drs. S. ©. — Maracaibo, Vene- 
zuela c/o Compania Shell de Venezuela 
Ltd, Apartodo 19 (g) 

COHEN, A. — Dwetta, W. Pakistan, c/o Photo- 
graphic Survey Corp. Litd., P. O. Box 76 (g8) 

COSTER, Dr. H.P. — Bellaire, Texas, 144 Whipple 
Drive (g) (gk) 

GROEN, geol. drs. H. A. — Kalulushi (Kitwe), 
Northern Rhodesia, 7 Davisstreet (8) 
GRIJZEL, P. VAN — Alphen a/d Rijn, Oranje- 

laan 39, (bg) (L.G.V.) 
TACOBSFGL I Zutphen, Weg naar Laren 13 


(bg) (U.G.V.) 

KNAAP, W. A. — Amstelveen, Troelstralaan 34 
(bg) 

OCRELORN Gr 2 Delft, van Leeuwenhoek- 


singel 28 (bg) (M.V.D)) 

WEGEN, G. VAN DER — Tandjong Pandan, 
Billiton, Indonesia, p/a G.M.B.-Billiton 
ZEE, gel. drs. T. J. VAN DER — South Porcupine, 

Ontario Canada, Powell Street 59 (g) (gk) 
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L. Bardowicks, 
Pomonoplein 67, Den Haag . 


STRECKENAUSBAU AUS VBS-STAHL 


(Vergüteter Bergbau-Sonderstahl) 


Im Jahre 1951 wurde im neu erbauten Werk 
der Gerlach G.m.b.H., Gladbeck i.W., mit der 
Vergütung von Thomas-Stahl für Strecken- 
ausbau (VBS-Stahl) begonnen. Die Festig- 
keitswerte betragen nach der Vergütung (Mit- 
telwerte aus über 300 Schmelzen): 


Bruchfestigkeit 75,— kg/mm? 
Streckgrenze 55,— 2 
Dehnung 20 %. 


Durch die Erhöhung der Festigkeitswerte 
ist es möglich, mit leichteren Profilen, als 
üblich, auszukommen. 


Die Vergütung verbessert ausserdem die 
Zähigkeit des Materials ganz ausserordentlich. 
Dies macht sich in den hohen Kerbschlag- 
zähigkeitswerten von über 5,0 mkg/cm? im 
gealterten Zustand bemerkbar. Bögen aus 
VBS-Stahl werden daher auch bereits seit 
mehreren Jahren in grossem Umfange kalt 
gerichtet. 


Bei der Fertigung duchlaufen die auf pas- 
sende Länge geschnittenen und gelochten 
Walzstäbe einen selbsttätig arbeitenden Hub- 
balkenofen, an dessen Ende sie auf eine 
Temperatur von etwa 900° C gebracht werden. 
Über einen Rollgang gelangen sie anschlies- 
send zur Biegevorrichtung, in der sie in 
einem Arbeitsgang gebogen und durch 
Absenken der gesamten Biegevorrichtung im 
Tauchbad abgeschreckt werden. Die Bogen- 
stücke bleiben während des Abschreckens in 
der geschlossenen Matrize, so dass sie sich 
nicht verwerfen können. 


Die Bogenteile werden nach dem Abschrek- 
ken zum Anlassofen gebracht, der ebenfalls 
als Hubbalkenofen ausgebildet ist. Die Tempe- 
raturen des Anlassofens richten sich nach den 
gewünschten Festigkeitswerten und der Zu- 
sammensetzung des Materials. Nach dem 
Austritt aus dem Anlassofen laufen die Bogen- 
stücke über ein Kühlbett, wo sie an der Luft 


abgekühlt werden. Zur Beschleunigung der 
Kühlung kann zusätzlich durch einen Ven- 
tilator Luft aus einem Röhrensystem auf die 
Bogenteile geblasen werden. 

Beide Öfen sind mit Temperaturreglern ver- 
sehen, die die gewünschten Temperaturen 
innehalten. Die Temperaturen werden laufend 
von Überwachungsgeräten aufgeschrieben, so 
" dass jederzeit nachgeprüft werden kann, ob 
\ die Bogenteile die richtige Wärmebehandlung 
erfahren haben, 
| Alle Bogenteile erhalten zur Sicherung 
“ gegen Herausfallen der Holzbolzen nach dem 
ICYC-ARC-Verfahren angeschweisste Stifte, 
von denen zurzeit monatlich rd. 50’000 ver- 
schweisst werden. 

Der Ausbau wird in starrer Ausführung 
mit Klammerlaschen und Schrauben, in ge- 
lenkiger Form mit Rollkeilgelenken und Holz- 
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boormachine 


voor ononderbroken 


läuferschalen und schliesslich nachgiebig mit 
Gleitfüssen geliefert. 

Bei der nachgiebigen Ausführung werden 
Klemmbacken, die einen Einsinkwiderstand 
von rd. 15 t erzeugen, bereits in der Werk- 
statt angeschraubt. Die Klemmkraft wird 
durch Eindrücken gehärteter Schneiden in den 
Steg der Bogenprofile erzeugt, so dass sie 
fast unabhängig von der Grösse des Anpress- 
druckes ist und jederzeit — auch nach längerer 
Standdauer und unabhängig von dem Grad 


‚des Rostens — anspricht. Die Klemmverbin- 


dung kann in der Grube verstellt, durch Lösen 


der Schrauben leicht geraubt und immer 


wieder ohne zusätzliche Teile wiederverwen- 
det werden. 

Das untenstehende Bild zeigt eine Anwen- 
dung der Gleitfüsse bei Türstockausbauten in 
Abbaustrecken der flachen Lagerung. 
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HEERLEN 


Bouwt sinds meer dan 60 jaar uitsluitend mijnbouwmachine 


GROSSLOCHBOHREN. 4%, 
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Großloch-Bohr- 
maschine „GB10" 
mit Druckluft- oder 
Elektroantrieb. 

Kontinuierlicher 
Vorschub = wirt- 
schaftliches Bohren 
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Vertegenwoordiging voor Nederland P. Bergstein, Handelmaatschappij C. V. 
Poststraat 19, Kerkrade. 


MIJNLOCOMOTIEF 


met .een 6-cylinder 
KROMHOUT 
DIESELMOTOR 
van 65 epk. 
in gebruik bij 
de Nederlandse 
mijnen 


%* 


AMSTERDAM 
Postbus: 959 Tel: 61611 


wessen LÖBBE-SCHAAFPLOEG 


DE NIEUWE KOOLWINNINGS-, LAAD- EN 

TRANSPORT-MACHINE. MET EEN CAPA- 

CITEIT TOT 1000 TON PER DIENST VOOR 

LAGEN VAN 40 CM EN HOGER, BIJ EEN 

HELLING VAN 30° BENEDENWAARTS 
EN 16° OPWAARTS. 


52/& H. v. Neindor! 


GEWERKSCHAFT EISENHÜTTE 


WESTFALIA LUNEN 


VERT. INGENIEURS-BUREAU „FERRUM” HEERLEN 


MASCHINENFABRIK GLÜCKAUF - GELSENKIRCHEN 


Dubbelwerkende Aluminium Omdrukcylinders voor 
kettingtransporteurs 


1 
E 30 mm P bij AN 530 kg \ 400 mm 


200° mm 21300 kg, ladlsuge pe 


gewicht 130/400 19 kg - hoogte 155 mm - breedte 230 mm 
gewicht 130/800 28 kg - hoogte 155 mm - breedte 230 mm 
gewicht 200/400 32 kg - hoogte 260 mm - breedte 250 mm 
gewicht 200/800 42 kg - hoogte 260 mm - breedte 250 


Uitschuifbare aluminium schoorstempels 


Vertegenwoordiging: 


INGENIEURSBUREAU „LIMAHA” U 


KERKRADE - ABTENLAAN 58 
TEL. K 4445 - 2156 


BOCHUM 
Vert. Ing. Bur. „FERRUM“, Heerlen 
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CLETON 


NIEUWE HAVEN 143-147. SCHIEDAM 


voor warmte- en koudeverliezen, 
TEL. 65111 (3 LIINEN) 


geluiddemping en trillingen. 
onbrandbaar, onaantastbaar voor 
N.V. Internationale Technische Handelmij 


vocht en schimmel. 
British Betting & Asbestos Ltd 


Alleenvertegenwoordigers van de 
Cleckheaton, England 
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transportbanden, op p.v.c. basis 
Telet 40767 
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H.R. SMITH N.V. 


KEIZERSGRACHT 520 - A'DAM - TEL.42012-41801 


Metaalgelijkrichtereenheden 


Agente en verlegenwoordigsier 


WESTINGHOUSE 


BRAKE AND SIGNAL COMP, LTD, 


steeds 
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dienst 


voor alle doeleinden 


Licentie voor Nederland, Indonesi& en Overzeese Rijksdelen 


Geschikt voor continu bedrijf 


GAPACITEIT cn 7 am. vegenerux 


9 Um 90 m3/min., vrii afgegeven hoeveelheid lucht 


KRACHTVERBRUIK 


0,0797 .0,0827 kWh. per m3 vrij afgegeven hoeveelheid lucht 


GOMPAGTE BOUW 
PERFEGTE UITBALANGERING 


UITSTEKENDE REFERENTIES 
KORTE LEVERTIIDEN 


type | m’/min, | type | m’/min, 
AR 1 93 ars | 283 
R3| 162 JaRz | 50 
RA, 215 Jana | 9 


Vraagt Allas om 


technische voorlichting op het 
gebied van luchtcompressoren 

luchtdrukwerktuigen 
en verfspuitinstallaties 
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BREEVAARTSTRAAT 48 - ROTTERDAM _ TEL 35191 
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W. A. HOEK’s 


MACHINE- EN ZUURSTOFFABRIEK 
Postbus 78 - Schiedam 


Ten 


5 


SIEMENS 


ELECTRICITEIT IN DE MIJNEN 


Centrale bediening De Siemens-bedieningslessenaar bevat centraal alle voor de ophaal- 
_ machine vereiste bedienings- en contröle-inrichtingen, zoals: 


1 bedieningshefboom 

2 hefboom voor de bedrijfsrem 

3 hefboom voor de noodrem 

‘4 hefboom voor personenvervoer 
5 diepteaanwijzer 

6 snelheidsregelaar 

7 contröle-instrumenten voor de rem en de ophaalmotor 


8 snelheidsmeter 


9 aanwijsinstrumenten voor de schachtsignaalinstallatie 
10 luidspreker en microfoon 


Alle instrumenten zijn hierdoor zeer overzichtelijk voor de ophaalmachinist 
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en maken een continu-contröle van de bediening mogelijk. 
NEDERLANDSCHE SIEMENS MAATSCHAPPI) N.V.- RUNSTRAAT 24 - s’GRAVENHAGE: TEL. 723810 
ALLEENVERTEGENWOORDIGING VAN! 
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S. K. A. ventilatormotoren 6000 pk STAATSMIJNEN Omzetters ZOUTINDUSTRIE 


Voor bedrijfszekerheid HEEMAF 


